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ABSTRACT 


Variable-Domain Displacement Transfer Functions were formulated for shape predictions of complex 
wing structures, for which surface strain-sensing stations must be properly distributed to avoid jointed 
junctures, and must be increased in the high strain gradient region. Each embedded beam (depth-wise cross 
section of structure along a surface strain-sensing line) was discretized into small variable domains. Thus, 
the surface strain distribution can be described with a piecewise linear or a piecewise nonlinear function. 
Through discretization, the embedded beam curvature equation can be piece-wisely integrated to obtain the 
Variable-Domain Displacement Transfer Functions (for each embedded beam), which are expressed in 
terms of geometrical parameters of the embedded beam and the surface strains along the strain-sensing line. 
By inputting the surface strain data into the Displacement Transfer Functions, slopes and deflections along 
each embedded beam can be calculated for mapping out overall structural deformed shapes. A long tapered 
cantilever tubular beam was chosen for shape prediction analysis. The input surface strains were 
analytically generated from finite-element analysis. The shape prediction accuracies of the Variable- 
Domain Displacement Transfer Functions were then determined in light of the finite-element generated 
slopes and deflections, and were found to be comparable to the accuracies of the constant-domain 
Displacement Transfer Functions. 


NOMENCLATURE 

a coefficient of first term of quadratic function, s(^) = a + b^ + , in./in. 

b coefficient of second term of quadratic function, £{<^) = a + b^ + , 1/in. 

c depth factor (vertical distance from neutral axis to lower surface of uniform embedded beam), 

in. 

c(x) depth factor (vertical distance from neutral axis to lower surface of nonuniform embedded 

beam at x), in. 

Cj = c(Xj), lower depth factor (vertical distance from neutral axis to lower surface of 

nonuniform embedded beam at x = x-), in. 

Cl upper depth factor (vertical distance from neutral axis to upper surface of nonuniform 

embedded beam at x = x, ), in. 

Cq value of c, at fixed end (embedded beam root), x = Xg = 0 , in. 

c„ value of c, at free end (embedded beam tip), x - x„ - I , in. 

DFF design limit load 

e coefficient of third term of quadratic function, £{^) = a + b^ + , 1/in^ 

hi depth of embedded beam at x = X, . 

i = 0,1,2,3,....,« , strain-sensing station identification number 

j dummy index 

k dummy index 

I length of embedded beam, in. 

n index for the last span-wise strain-sensing station (or number of strain- sensing domains) 

R radius of curvature, in. 

RRF rotated reference frame, aligned with wing’s out-of-plane direction 

SPAR Structural Performance And Resizing 

X, y Cartesian coordinates (x in embedded beam axial direction, y in lateral direction), in. 

X, axial coordinate (or symbol) associated with /-th strain-sensing station, in. 


1 


y(x) embedded beam deflection at axial location x, in. 

= j(X;) , embedded beam deflection at axial location, x - x^, in. 

(Ji)pred. predicted value of J, , in. 

(Ji )sPAR value of calculated from SPAR program, in. 

AZ = (x, - x,_j )-lln, constant-domain length (distance between two adjacent strain-sensing 

stations {X;_j , x, } , in. 

(A/), = (X; -x,_i) , variable-domain length (distance between two adjacent strain-sensing stations 

{x,._i,x,.}, in. 

£{x) lower surface bending strain at axial location X, in/in 

£. = £(x.), lower surface bending strain at strain- sensing station, x, , in/in 

£. upper surface bending strain at strain-sensing cross section, x, , in/in 

0(x) embedded-beam slope in reference to x-system, deg 

6*, = 6*(x, ) , slope at strain- sensing station, x - x^, deg 

(^i)pred. predicted values of 6*; , deg 

(6*, ) SPAR SPAR-calculated values of 6*, , deg 

^ = X - x,_[ , local axial coordinate measured from strain-sensing station, x,_j , in. 

INTRODUCTION 

For structural deformed shape predictions using distributed surface strains, the Displacement Transfer 
Functions (refs. 1-9) are needed to convert the surface strains into out-of-plane deflections for mapping out 
the overall structural deformed shapes. The surface strain data are to be obtained at the strain- sensing 
stations discretely distributed along the strain- sensing line lying on the surfaces of a planer structure (for 
example, an aircraft wing). The depth-wise cross section of the structure along the strain- sensing line can 
be considered as an embedded beam with width of the strain-sensing line (fig. 1). 

In the formulations of earlier Displacement Transfer Functions (refs. 1-9), the embedded beam was 
first discretized evenly into multiple small domains (constant-domain lengths) with strain- sensing stations 
located at the domain junctures. Thus, within each small domain, the depth variation of the embedded beam 
could be described with a linear function, and the surface strain distribution could be described with either 
linear or nonlinear function. Such a piecewise approach enabled piecewise integrations of the embedded- 
beam curvature equation (second order differential equation) to yield beam slope and deflection equations 
in recursive forms. Each set of recursive slope and deflection equations were then combined into a 
summation- form deflection equation (called the Displacement Transfer Function) for each embedded beam. 
In the Displacement Transfer Function, the deflections are expressed explicitly in terms of geometrical 
parameters of the discretized embedded-beam and surface strains along the strain-sensing line. By inputting 
the surface strains into the Displacement Transfer Function, one can then calculate deflections along the 
embedded-beam for generating the elastic curve. By using multiple strain-sensing lines, one can then map 
out all the elastic curves of the multiple embedded beams to construct the overall structural deformed shapes 
(under bending and torsion) for visual display. The Displacement Transfer Functions combined with a 
strain- sensing system, thus created a revolutionary new structure-shape-sensing technology, “Method for 
Real-Time Structure Shape-Sensing, ” (U.S. Patent Number 7,520,176, ref. 6), which is very attractive for 
application to the in-flight deformed shape monitoring of flexible wings and tails of unmanned flight 
vehicles by the ground-based pilot for maintaining safe flights. In addition, the real-time wing shape 
monitored could then be input to the aircraft control system for aero-elastic wing shape control. 
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The shape prediction accuracies of the earlier constant-domain Displacement Transfer Functions (based 
on piecewise-linear strain representations) (refs. 1-7) were successfully validated analytically using finite- 
element analyses of different sample structures such as cantilever tubular beams (uniform, tapered, slightly 
tapered, step- wisely tapered), two-point supported tapered tubular beams, flat panels, and tapered wing 
boxes (un-swept and swept). 

Also, by using piecewise-nonlinear strain representations in the formulations (ref. 8), the shape 
prediction accuracies could be improved greatly, especially for highly tapered beam structures with rapidly 
changing strain gradients. The Displacement Transfer Functions were originally formulated for the straight 
embedded beams. However, with the introduction of empirically established curvature-effect correction 
terms, those Displacement Transfer Functions could be used for shape predictions of slender curved 
structures with different arc-angles up to 360 degrees, a complete circle (ref. 10). All the earlier 
Displacement Transfer Functions (refs. 1-9) were formulated based on uniform discretization of the 
embedded beam into multiple domains with equal domain lengths. 

For the shape predictions of certain long complex wing structures consisting of multiple sections joined 
together, one must avoid locating the strain-sensing stations right at the jointed junctures, where the strain 
outputs could be erroneous. Also in the region where the strain gradient changes rapidly, shorter strain- 
sensing intervals (shorter domains) are needed to improve the shape prediction accuracies. Such situations 
demand the use of variable-domain lengths for properly locating the strain-sensing stations. Therefore, 
variable-domain Displacement Transfer Functions are needed for deformed shape calculations of a complex 
structure for which a non-uniform discretization approach is required. 

This technical publication presents the formulations of the variable-domain Displacement Transfer 
Functions (unevenly distributed strain-sensing stations) by extending the previously developed constant- 
domain Displacement Transfer Functions (evenly distributed strain- sensing stations) (refs. 1-9) for 
structural deformed shape predictions. 

A long tapered cantilever tubular beam was chosen for the shape prediction accuracy studies of the 
variable-domain Displacement Transfer Functions. Finite-element analysis was used to analytically 
calculate the input surface strains, and to generate reference yardsticks (slopes and deflections) for the shape 
prediction error analysis. The variable-domain Displacement Transfer Functions were found to provide 
highly accurate shape predictions, similar to the constant-domain Displacement Transfer Functions. 

NONUNIFORM AND UNIFORM BEAMS 

In the present technical publication, the term “beam” implies the embedded beam, which is defined as 
the depth-wise cross section of the structure along the strain- sensing line, and different from the isolated 
classical Euler-Bernoulli beam. The nonuniform beam implies the embedded beam with varying depth. The 
slightly tapered beam implies the embedded beam with slowly changing depth. The uniform beam implies 
the embedded beam with constant depth. 

BASICS OF CONSTANT-DOMAIN DISPLACEMENT TRANSFER 

FUNCTIONS 

This section reviews the basics of formulations of the earlier Displacement Transfer Functions for 
constant-domain lengths (refs. 1-9). 
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Curvature- Strain Equation 


Formulations of the earlier constant-domain Displacement Transfer Functions for structure deformed 
shape predictions stemmed from the integrations of the following geometrically established shifted 
curvature-strain differential equation for the embedded beam (ref. 9): 

dx^ c(x) R 

In equation (1), x is the undeformed axial coordinate, y is the lateral deflection, s(x) is bottom-surface 
bending strain, and c(x) is the embedded beam depth factor, which is defined as the vertical distance from 
the neutral axis to the bottom surface of the embedded beam. Equation (1) is the simplified Lagrangian 
curvature equation, 1/ R(x) = (d^y / dx^)j -sjl-idy / dx)^ , with the term, (dyldx)^ , neglected. Neglection 

of the (dy/dx)^ term is equivalent to setting axial displacements to zero (called shifted Lagrangian 
formulation) (refs. 9, 11, and 12). Keep in mind that equation (1) is referred to as the undeformed 
x-coordinate, and is not the simplified form of the classical curvature equation, 

2 2 / 2 3/2 

llR = (d y I dx )/[l + {dy / dx) ] , which is referred to as the deformed x-coordinate (Eulerian 

formulation) (refs. 9, 11, and 12). 

The seven types of constant-domain Displacement Transfer Eunctions formulated earlier, based on 
equation (1), were found to be quite accurate [0.0235~0.3852 percent error range at beam tip] in the shape 
predictions of a tapered cantilever tubular beam under both small and large bending deformations, even 
with the free-end deflection reaching as large as 94 percent of the structural length (ref. 9). 


Discretizations 

To integrate equation (1), the embedded beam of length, I, lying along the strain-sensing line (fig. 1), 

was first discretized evenly into n number of small domains, x,_j < x< X- (/ = 1,2,3, ,n), with constant- 

domain lengths, M = l ! n[= (x- -x,_j)], so that within each small domain, x,_j <X<X-, the beam depth 
factor, c(x) (vertical distance from neutral axis to lower surface of beam at axial location, x), could be 
described with linear functions, and the surface strain, s(x ) , could be described with either linear or 
nonlinear function as shown in equations (2)-(4a): 


Piecewise-linear depth factor representations (refs. 1-7): 


c(x) = c,._i-t(c,.-c,_i)- 


M 


; (x,._i <x<x,.) 


Piecewise-linear strain representations (refs. 1-7): 


S{x)=€i_^+(S,-S,^^) 


X-X. 


i-l 


A/ 


; (x,._i <x<x,.) 


( 2 ) 


( 3 ) 
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Piecewise-nonlinear strain representations (ref. 8): 




2 A/ 


2(A/r 


(4a) 


(x,._i <x<x,.) 


Note that, when i = n, equation (4a) eontains a nonexistenee (fictitious) strain, , which can he 
estimated hy using the following three -points extrapolation equation: 

^,,+1 = ^n-i - 3^„-i + 3^„ (4h) 

Equation (4h) was obtained from equation (4a) hy setting i = n - \ and setting the term (x - x,_j ) to 
(^«+i - ^«-i) = 3 A/ in the region (x„_i < x < x„^i). 

The piecewise approach enabled the piecewise integrations of the curvature-strain differential equation 
(1) to yield closed-form beam slope and deflection equations in recursive forms. The recursive slope and 
deflection equations were then combined into one deflection equation in dual summation forms (called the 
constant-domain Displacement Transfer Function), which are expressed in terms of embedded beam 
geometrical parameters and surface strains obtained at evenly spaced strain- sensing stations (at domain 
junctures) along the strain- sensing line. By inputting the surface strain data into the Displacement Transfer 
Function, one can then calculate the out-of-plane deflections along each embedded beam. By using multiple 
strain sensing lines, one can then map out the overall structural deformed shapes (under bending and 
torsion) for visual display. 


LISTS OF CONSTANT-DOMAIN DISPLACEMENT TRANSFER 

FUNCTIONS 

This section lists previously formulated key Displacement Transfer Functions for constant-domain 
lengths, M = Hn (uniform discretization) (refs. 1-9). Those constant-domain Displacement Transfer 
Functions can then be used as a basis to formulate the variable-domain lengths Displacement Transfer 
Functions (nonuniform discretization) presented in the subsequent sections. Depending on the structural 
types, one can choose a suitable Displacement Transfer Function from the following list for shape 
predictions. 


Based on Piecewise-Linear Strain Representations 

The following two sets of slope and deflection equations were formulated using piecewise-linear 
functions [eqs. (2) and (3)] to describe the distributions of both depth factor and surface strain, 
{c{x),£{x)]. 

1. Nonuniform Embedded Beams (Linear Strain) (Ci ^ c, ) (refs. 1, 3) 

Slope equation: 


tan 9: - A/ 


r - £. 


Cm -c,' 




(Cm -c, ) 


Ci 


-t tan 


1-1 


(i= 1,2,3, ,n) (5a) 
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Deflection equations: 
a. In recursive form: 


= {Mf 




2(c,.-i-c,.) (C;^i-C,.) 


c, log— ^ + (c,.^i -c,.) 
C,- , 


■ + + Al tan (5b) 

(/' = 1,2,3, ,n) 


b. In dual summation form: 




M 2(c^._i-cp 




■^jr 




c.log, + (c, i-c,) 




t-1 






Contributions from deflection terms 


fill fl + fiiifi flf^log^ 








■+ Jo +(OA/tan6*o 

' " V ^ 

=0 for cantilever beams 


Contribution from slope terms 


U - 1,2,3, ,n) 


(5c) 


Equations (5) are called tbe Constant-Domain Nonuniform Displacement Transfer Functions for a 
nonuniform embedded beam based on piecewise-linear strain representations. The term “Displacement 
Transfer Function” was used because equation (5c) transforms the surface strains into out-of-plane 
deflections for plotting the deformed shape of the embedded beam. 


Note that equations (5) are not applicable to the limit case of uniform beam (C; j = c, = c) because of 
mathematical indeterminacy (0/0) caused by log(Cj/C(_i) ^ 0 and (c,. - c,._j) ^ 0 . However, by expanding 
the logarithmic terms in equations (5) in series forms in the neighborhood of (c-_^ /c.)— > 1, the factor 
(c, — C._[) in the denominators can be eliminated (refs. 1, 3). The resulting Fog-expanded equations (6) 
listed below are for slightly nonuniform embedded beams [(c,_j / c,) — > 1], and applicable to the limit case 
of uniform beams (c, j = c, = c) . 


2. Slightly Nonuniform Embedded Beams (Linear Strain) [(c, j /c, ) — > 1] (refs. 1, 3) 
Fog-expanded slope equation: 


tan 0. - 


Al 


‘ 2c, , 




-l-tan6';_i ; (1=1,2, 3, ,n) 


(6a) 


Fog-expanded deflection equations: 
a. In recursive form: 


3'i = 


(Aiy 

6c,. 


[f \ 


- C; 


3 !- 

^i-1 + 

_V ^i-1 y 



-t y,^i -r AZ tan ; (Z=I,2,3, ,n) (6b) 
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b. In dual summation form: 


>'/=■ 


{Mf 


6 


t- 




A ^ 

3--^ 

V y 




” + >0 


Contribution from deflection terms 


+ - 


x2 i-l 


,/=l 


f.- M 

■/ X 


1 

(i-j) 

2--fM 

£ ,+£■ 


c , 

c , 

J 1 J 

1 

L 2-1 [ 

. V j-i y 


J 


> + (i)A/tan0o 


Contributions from slope terms 


(6c) 


Equations (6) are called the Constant-Domain Log-expanded Displacement Transfer Functions for 
slightly nonlinear (including uniform) embedded beams based on piecewise-linear strain representations. 


Based on Piecewise-Nonlinear Strain Representations 

The following two sets of slope and deflection equations were formulated by using the piecewise-linear 
function [eq. (2)] to describe the actual distribution of the depth factor, c(x), and the piecewise-nonlinear 
function [eq. (4)] to describe the actual distribution of the surface strain, s(x ) . 


1. Nonuniform Beams (Nonlinear Strain) (C; j ^ c, ) (ref. 8) 


Improved slope equation: 

^ A/ 

tan 0: - 


2(C -c,._i)' 

A/ 

'4(c,-c,_i)^ 


- [(2C; - c,._i )(c,.f ,._1 - ]log, 

[(5c,. -3 c,._i)£-,._i -2(3c,. -c,._i)£-,. -t(c,. -tc,._i)£-,._^i]-ttan^,._ 


(7a) 


(i = 1,2,3,...., n) 


Improved deflection equations: 
a. In recursive form: 


J/ = 


{Mf 

2(c,.-c,_,y 

(Mf 

12(c,.-c,._i)' 


-[(2c, -c,._i)(c,.^,._i -2c,._if,.) + c,c,._if.,i] 


c,log,^-(c,.-c,._i) 


(7b) 


-[(8c,. - 5c;_i)f,._i - 2(5c,. - 2c,._i)f,. + (2c,. + c,_i )f,.,, ] + y,._i + A/ tan 

(i = 1,2,3,...., n) 
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b. In dual summation form: 





c, 


X 




Contribution from deflection terms 



(7c) 


+ (/)A/ tanSg 


Contributions from slope terms 


{i - 1,2,3, ,n) 


Equations (7) are called the Constant-Domain Improved Displacement Transfer Functions for nonuniform 
embedded beams based on piecewise-nonlinear strain representations. 

Equations (7) are not applicable to the limit case of uniform beam ( c,_j - - c ) because of 

mathematical indeterminacy (that is, 0/0). To handle the limit case of uniform beams (c, j = c, = c), the 
logarithmic terms in equations (7) could be expanded in series forms in the neighborhood of (c,_j / c, ) ^ 1 
(ref. 8), and obtain the following Eog-expanded equations (8a), (8b), and (8c) for slightly nonuniform 
beams, and which are applicable to the limit case of uniform embedded beams (C; j = c, = c) . 

2. Slightly Nonuniform Embedded Beams (Nonlinear Strain) [(c, j / c, ) — > 1] (ref. 8) 

Eog-expanded slope equation: 



+ tw9^_^ (8a) 


(/ = 1,2,3,. ...,n) 


Eog-expanded deflection equations: 
a. In recursive form: 



(8b) 


+ y,._j -I- A/ tan 


(/ = 1,2,3,. ...,n) 


b. In dual summation form: 





(8c) 


Contributions from deflection terms 



+ (/)A/ tan0(, 


Contributions Irom slope terms 


(z = 1,2,3, 


Equations (8) are called the Constant-Domain Log-expanded Displacement Transfer Functions for 
slightly nonuniform (including uniform) embedded beams based on piecewise-nonlinear strain 
representations. 


This section describes the formulations of the Displacement Transfer Functions for variable-domain 
lengths (nonuniform discretization cases). Let the embedded beam be discretized into n domains with 
variable-domain length, (A/);(=x, -x,_j) (z = 1,2,3,. ...,zz) [A/ in figure 1 replaced with (A/)J, then the 
beam depth factor, c(x) , within each small domain, X;_j < x < X; , between the two adjacent strain-sensing 
stations, {x,_j ,xj, can be represented with the piecewise-linear function, and the surface strain, s(x) , 
represented with either the piecewise-linear or the piecewise-nonlinear function (quadratic function). 

Piecewise-linear depth factor representation as shown in equation (9): 


FORMULATION OF VARIABLE-DOMAIN DISPLACEMENT 

TRANSFER FUNCTIONS 


X — X ■ 

c(x) = C;_1 -I- (c,. - C,._1 ) ; (x,._i < X < X,. ) 


(A/),. 


(9) 


Piecewise-linear strain representation as shown in equation (10): 


.X . x-x 

£(X) = ZT,._i -I- (£i - £i_i ) — ^ 


(x,._i <X<X,.) 


(10) 
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Piece wise-nonlinear strain representations (fig. 2) (see Appendix A for derivation): 


£(x) = 



(11) 


(x,._i <x<x,.) 


When i = n, equation (11) contains a nonexistence (fictitious) strain, By setting i = n-l and 
setting the term (x-x,_i) to (x„_,i - x„_i) = [(A/)„_i -t (A/)„ -t (A/)„_^J in the region (x„_i < x < x,,^;) , 
equation (11) becomes the following three-points extrapolation formula for the calculation of the fictitious 
strain, 


For a special case of constant-domain lengths, (A/)j = (A /)2 = (A /)3 = - AI , it is easy to show that 

equations [(1 1), and (12)] will degenerate respectively into equations [(4a), and (4h)] (see Appendix A). 


The variahle-domain slope and deflection equations in recursive forms can he obtained from the 
constant-domain recursion slope and deflection equations [(5a), and (5b)-(8a), and (8b)] by simply 
changing the constant-domain length, A/(=//n), into variable-domain lengths, (A/),(= X; -x,_i) . 
However, the derivations of the variable-domain Displacement Transfer Functions (recursive slope and 
deflection equations combined into dual summation forms) require considerable mathematical 
manipulations (see Appendix B). 


The following two sets of variable-domain slope and deflection equations were formulated by using 



( 12 ) 


LISTS OF VARIABLE-DOMAIN DISPLACEMENT TRANSFER 

FUNCTIONS 


Based on Piecewise-Linear Strain Representations 


piecewise-linear functions to describe the piecewise variations of the depth factor and surface strain, [c(x) 
s{x)}, [eqs. (9), and (10)]. 
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1. Nonuniform Beams (Linear Strain) [(A/), ; ^ ( A/); ; T^c,.)] 

Slope equation [see eq. (5a)]: 


tan^,. = (A/),. 


^/-l I ^i-Ni 

Ci-l-Ci 




log— 


+ tan 0- 


/-I 


(13a) 


(i = 1,2,3, ...„n) 


Deflection equations: 

a. In recursive form [see eq. (5b)]: 




Si-Ni-SNi-i 


2 (cm-c,) (c,.^i-c,.)^ 


c, log^ + (c,._i -c,.) 


■'i-l 


+ J/- 1 +(A/), tan 6 ',._i (13b) 

(i = 1,2,3,. ...,n) 


b. In dual summation form (see Appendix B for derivation): 


,;=i 


£; ,c. -e.c. , / c. ' 

^ — + ("m - S' ) 


' + )'o 


Contributions from deflection terms 


y=i 


X (A/), 


k=i+\ 


(A/),. 


^j-i , ^j-Nj ^fj-i S 


t-i-S (S-i-S) 


'y-i 


I (AO, 


;=i 


tanSn 


Contributions from slope terms 


(13c) 


(z = 1,2,3,... .,zr) 


Equations (13) are called the Variable-Domain Nonuniform Displacement Transfer Functions for 
nonuniform embedded beams based on piecewise-linear strain representations. For constant-domain 

lengths, (A/)j = (A /) 2 - (A /) 3 = = A/ , equation (13c) degenerates into equation (5c). Equations (13) are 

not applicable to the limit case of uniform beam (c, j = c, = c) because of mathematical indeterminacy (that 
is, 0/0). However, by expanding the logarithmic terms in equations (13) in series forms in the 
neighborhood of (c, j /C/)—>l (refs. 1, and 3), one can obtain the following Fog-expanded equations (14) 
for slightly nonuniform embedded beams [(c, j / c, ) — > 1 ], which are applicable to the limit case of uniform 
beams (c, , = c, = c) . 


2. Slightly Nonuniform Embedded Beams (Linear Strain) [(A/); j ^(A/),;(c,/c,_j)^l] 
Slope equations: 


tan 


(A/), 


f ^ 


C: 


2 !- 


_V y 



+ tan 6^_^ 


(i = 1,2,3,. ...,zz) 


(14a) 
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Deflection equations: 
a. In recursive form: 




(A/)f 


6c, -1 

b. In dual summation form: 


e-tJ' 


+ y,._i + A/ tan6*,._i ; (i ~ 1,2,3,.— ,«) 


(14b) 


6 C ^-1 


/ A 

3--^ 

V ; 


e._,+8j 


■ + Xo 


Contributions from deflection terms 


+ - 


tS- 


2-tr 


t (AZ)* 

(Al), 




> 

1 

_k=J+] 

^j-i 

A A-i J 



_ t=i 


tan0„ 


Contributions from slope terms 


(14c) 


(i - 1,2,3, ,n) 


Equations (14) are the Variable-Domain Log-expanded Displacement Transfer Function for slightly 
nonuniform (including uniform) embedded beams based on piecewise-linear strain representations. The 
derivation of equation (14c) is similar to that of equation (13c) (see Appendix B). For constant-domain 

lengths, (A/)j = (A /)2 = (A /)3 = - Al , equation (14c) degenerates into equation (6c). It must be 

mentioned that the shape prediction accuracy of equation (14a), (14b), and (14c) was experimentally 
validated (see section: Experimental Validations of Shape-Prediction Accuracies). 


Based on Piecewise-Nonlinear Strain Representations 

The following two sets of slope and deflection equations were formulated by using the piecewise- 
nonlinear function to describe the variation of the surface strain, £{x) [eq. (11)], but the piecewise-linear 
function was used to describe the variation of the depth factor, c(x) [eq. (9)]. For simplifications, let the 
coefficients of the second and third terms of equation (1 1) be represented respectively with symbols A, and 
Bj as shown in equations (15) and (16) respectively: 


A, =- 


(A/),,i[2(A/), +(AQ,,Jg,_i - [(A/), +(A/),,J^g, +(A/)fg,,, 
(A/),(A/),,i[(A/), +(A0,,J 


(15) 




(A/),,ig,_i - [(A/), +(A/),,Jg, -r(A/),g,. 
(A/),(AZ),,,[(A/), +(AZ),,J 


(16) 


Then, the improved slope and deflection equations can be written as (see Appendix C for derivations): 
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1. Nonuniform Embedded Beams (Nonlinear Strain) [(A/); j ^ (A/), ;(c,._i ^ C;)] 
Improved slope equation: 


tan 6*. = ■ 


(A/), 


+ 


(A/), 


- A,(AZ),c,_,(c, -c,_J + fi,(A/)fct Jlog^ 


A 

c,_ 


2(c.-c._J 


[2A,(A/),(c, -c,_,) + B,(A/)f(c, -3c,_i)] + tan6',_, 


(17a) 


(Z = l,2,3,....,n) 


Improved deflection equations: 
a. In recursive form: 


yi = 


(A/)? 


+ - 


(A/)f 

6(Ci-c,_iy 


- [(c,. - c,_i f - A, (A/),, (c, - c,._i )c,._, + B. (A / )f cf_^ ] 


c. 




(17b) 


- [ 3A( A/), (c, - c,_, ) + ( A/)f (c, - 4c,_, ), ] + y,_, + ( A/), tan 


(/ = 1,2,3, ...„n) 


b. In dual summation form: 


j=i 


{Alt 


{Cj-c,J 


[(c^ - A/A/).(c. + B,(A/)Jc^,]x 

c. 

{Cj-Cj-ft 


X 


'y-1 


+ 


(A/) 


6{cj-c._ft 


[3A,(A/),(c^. - c^_j) + fi,(A/)J(c^. - 4c^_,)] 




+I 

j=i 


i {AD, 


k=j+\ 


Contribution from deflection terms 
\2 , 


(A/). KCj Cy_j) £^._J Aj(Al)jCj_fCj Cj 

{cj-cj_,tyBfMtjcu 


j-i 


+- 


(AZ)^ 


2(c,-c,-ir 


2 A, ( AZ), (c, - c,_, ) + B. m] (c. - 3c ) 


Contributions from slope terms 


+ 


S(A0, 


J=1 


tan0„ 


Contributions from slope terms 


(17c) 


(Z = 1,2,3, ,n) 


Equations (17) are called the Variable-Domain Improved Displacement Transfer Functions for non- 
uniform embedded beams based on piecewise-nonlinear strain representations. Derivation of equation (17c) 
is similar to that of equation (13c) (see Appendix B). For the constant-domain lengths, 
(AZ)i = (AZ )2 = (AZ )3 = = AZ , equations (17) will degenerate into equations (7) (see Appendix C). 
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Equations (17) are not applicable to the limit case of uniform embedded beam (c,_j - - c) because 

of mathematical indeterminacy (that is, 0/0). Therefore, the logarithmic terms in equations (17) were 
expanded in series forms in the neighborhood of {c-_^l C-)^ \ , to yield the following Log-expanded 
equations (18) for slightly nonuniform embedded beams [(c,_j /c,) — > 1], applicable to the limit case of 
uniform beams (c, j = c, = c) (see derivations in Appendix D). 

2. Slightly Nonuniform Embedded Beams (Nonlinear Strain) [(A/); j ^ (A/), ; (c, / c,_,) — > 1] 
Log-expanded slope equation: 


tan 0 = 


(A/), (6cti + (c,. - c,._, )(2c. - 5c,._i )) 

l_-A,(A/),c,_,(2c, -5c,_J + 2B,(A/)fcti, 


-ttan6';_j (18a) 

(/ = 1,2,3,. ...,n) 


Log-expanded deflection equations: 
a. In recursive form: 


(A/)f 

12c, 


[6c, + c,._i(10c,. - 3c,._i)(c, - c,._i)- 3(c,^ - c,^i)]^,._, 
+A,(A/);C,._i [2c,^i + (3c,. - 4c,_i)(c; - c,._i)] 
+fi,(A/)fc,^,[c,_.-3(c,-c,_j] 


>-ty,._;-r(A/),.tan6',._i (18b) 

(/ = 1,2,3,. ...,«) 


b. In dual summation form: 


1 ^ m] 


12"^ r'* 

M ‘'i-l 


[6c'_i + c^._i(10c. - 3c,._i)(c. - Cj_,)-3(c] - 
-hA . (A/), [2c"_i -t c,._i (c . - c )(3c^. - 4c ,_i )] }■ -r 

+Z?,(A/)^c.?,[c^._,-3(c^,-c^._,)] 


Contributions from deflection terms 


(18c) 


1 


1 

< 

3 

(6c[-i + (Cj - c,_i)(2c. - 5c^_i ))£,._! 

+ 

< 

1 

tan0Q 

k= j+\ 


[-Aj c,_, (2c, - 5c,_, ) + 2B. (Alf. \ 


_ .'-1 



Contributions from slope terms 


(i = l,2,3,....,n) 


Equation (18) are called the Variable-Domain Log-expanded Displacement Transfer Lunctions for 
slightly nonuniform (including uniform) embedded beams based on piecewise-nonlinear strain 
representations. Derivation of equation (18c) is similar to that of equation (13c) (see Appendix B). Lor 

constant-domain lengths, (A/)j = (A /)2 = (A /)3 = = A/ , equation (18c) degenerates into equation (8c) 

(see Appendices D, and E). 
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CHARACTERISTICS OF DISPLACEMENT TRANSFER FUNCTIONS 


The Displacement Transfer Functions listed above are purely geometrical in nature, relating surface 
strains and depth factors to out-of-plane deflections, and contain no material/structural properties. Thus, in 
the shape predictions using the Displacement Transfer Functions, no knowledge of material/structural 
properties is required. In fact, the effects of material/structural properties are felt only by the surface strains 
and influence their outputs. 

Another key characteristic of the Displacement Transfer Functions is that the deflection, y- , at the 
strain- sensing station, x, , is expressed in terms of the inboard surface strains, at 

all the inboard strain- sensing stations, (Xg,XpX 2 ,....,X;), including the strain, Sj, at the current strain- 
sensing station, x^, where y, is calculated. The outboard surface strains, are not 

needed in the calculations of y, . 

DETERMINATIONS OF DEPTH FACTORS 

To use the Displacement Transfer Functions to convert surface strains into deflections, the depth 
factors. Cl (i = , of the embedded beam are needed in the deflection calculations. To determine 

the depth factors, c^ (that is, locating the neutral axis), one can use two strain- sensing lines, each of which 
are lying respectively on the lower and upper surfaces of the embedded beam. The lower and upper surface 
strains, {s-, €■} (i = l,2,3,....,n) , can then be used to determine the lower and upper depth factors, 

{c., C,.}. 

As shown in figure 3, R is the radius of curvature of the neutral axis of the embedded beam under 
upward bending; AqBq is the undeformed segment length along the neutral axis; CD is the elongation of 
the lower surface from un-deformed length, AB{= AqBq) to deformed length, AC; —CD is the contraction 
of the upper surface from undeformed length, AB{= AqBq) to deformed length, AC. Then, through 
geometrical similarities of the curved triangles, AO AC, ABqBC, ABqBC (fig. 3), the lower and upper surface 
strains, {s^, can be related respectively to the lower and upper depth factors, {c, , C;},as: 

BC BC c - -BC -BC c 

s. = = = — (tension) ; £■ = — = = — (compression) (19) 

‘ABA,B,R AB A,B, R^ ^ 

Keep in mind that in upward or downward bending, the signs of the lower and upper surface strains, 
{Sj, €■}, are always opposite. By eliminating R in equation (19), the two depth factors, {c, ,C^ } , can be 
related as: 





a = 1,2,3,....,«) (20) 
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The depth of the embedded beam, h., at strain-sensing station, jc,, is given by: 

h. = c,. -h Cl ; (i = 1,2,3, (21) 

In view of equations (20) and (21), one can establish the depth factor equations in the forms: 

Ci = — ; Ci = — ; (i = 1,2,3,. ...,n) (22) 

Equation (22) can also be used for downward bending, and is applicable to both solid beam cross sections, 
and hollow beam cross sections (for example, aircraft wing cross sections). 

If {£•,, Sj} have the same magnitudes (opposite signs), then equation (22) gives Ci = Ci = hi/2, 
indicating that the neutral axis is located at the beam half depth. Equation (22) was used with great success 
in the shape predictions of both Ikhana (General Atomics Aeronautical Systems Inc., San Diego, California) 
composite wings (66-ft wingspan) (ref. 3) and Global Observer (AeroVironment Inc., Monrovia, 
California) composite wings (175-ft wingspan) (ref. 13), for which the locations of neutral axes were 
unknown. Therefore, equation (22) had to be used to calculate the depth factors, {C;,g} ^ for each 
embedded beam using the associated surface strains, { . 

STRUCTURE USED FOR SHAPE PREDICTION ACCURACY STUDIES 

To compare the shape-prediction accuracies of the newly formulated variable-domain Displacement 
Transfer Eunction with the earlier constant-domain Displacement Transfer Eunctions, a sample structure 
was needed. The structure chosen for the prediction accuracy analysis was an aluminum tapered cantilever 
tubular beam with geometries listed in table 1 . 

Table 1. Geometries of aluminum tapered cantilever tubular beam. 


Z, in. 

t, in. 

Co , in. 

c„ , in. 

a{=tan ^[(cg -c„)/Z]}, deg. 

(Length) 

(Wall thickness) 

(Root depth factor) 

(Tip depth factor) 

(Taper angle) 

300 

0.02296 

4 

1 

0.5729 


The embedded beam (depth-wise structural cross section along the strain- sensing line) shown in figure 
1 can also represent the embedded beam of the tapered cantilever tubular beam. The embedded beam for 
the tubular beam is defined as the span-wise vertical cross section along the strain-sensing line which is 
coincidental with the bottom generatrix (a straight line for generating tubular surface through sidewise 

circular motion). Because the depth factor, Cj(t = 1,2,3, ,n) (fig. 1), is the known local radius, only one 

strain- sensing line is needed to obtain surface strains, £((t = 1,2,3, ,«) (fig. 1) for shape predictions. An 

upward point load of P = 300 lb (or P = 600 lb) was applied at the beam free end. 

In the present technical publication, sixteen (n - 16) strain-sensing domains (distance between two 
adjacent strain-sensing stations) were used for the constant-domain cases with equal domain lengths of 
AZ = Z/n = 300/16 = 18.75 in. To create variable-domain cases, the strain-sensing stations at i =4, 8, and 
12 of constant-domain cases were removed to double the domain length to (AZ); = 2 X 18.75 = 37.5 in. in 
those three regions. The rest of the domain lengths remain unchanged [that is, (AZ); = 18.75 in.] 
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ANALYTICAL SHAPE PREDICTIONS 


The shape prediction analysis in the present technical publication report is called analytical shape 
prediction analysis. Namely, instead of using experimentally measured surface strains, the Structural 
Performance And Resizing (SPAR) finite-element computer program (ref. 14) was used to analytically 
calculate the surface strains for input to the Displacement Transfer Functions for shape calculations. Also, 
for estimations of shape prediction errors of different Displacement Transfer Functions, the SPAR- 
generated slopes and deflections were used as reference yardsticks. 

FINITE ELEMENT ANALYSIS 

The purpose of the finite-element analysis is to generate input surface strains, and to calculate the slopes 
and deflections, which were used as the reference yardsticks in the shape prediction error analysis. 

Finite-Element Model 

The SPAR finite-element model generated for the tapered cantilever tubular beam has the following 
numbers of nodes and elements: 

3673 nodes = 3636 nodes for the tube wall (101 axial nodes x 36 circumferential nodes) -i- 
37 nodes on the beam-tip disk 

3636 four-node elements = 3600 elements for the tube wall + 

36 elements for beam-tip disk outer region 

36 three-node elements for beam-tip disk central region 

The upward point load of P = 300 lb (or P = 600 lb) was applied at the beam-tip disk central node of 
the SPAR model. 

Figure 4 shows the undeformed and deformed shapes of the SPAR model generated for the tapered 
cantilever tubular beam subjected to beam-tip load of P = 300 lb (orP = 600 lb). Note that for the loading 
case of P = 300 lb, the beam-tip deflection reached y„Kl41.6 in., and the beam-tip slope reached as large 
as 48 deg. For the doubled loading case of P = 600 lb, the beam-tip deflection was doubled to y„~ 
283.2 in. (because of linear elasticity), however, the beam-tip slope was not doubled, but only increased by 
1.375 times to 6*„ s; 66 deg (=1.375 x 48 deg). 

It must be mentioned that the SPAR model has 100 elements in the axial direction (fig. 3). For n = 16 
strain- sensing domains, the even number of strain sensing stations can be either coincidental with SPAR 
nodes or at the midpoints between two adjacent SPAR axial nodes. However, the odd numbers of strain- 
sensing stations are not precisely at the SPAR nodes. Therefore, the shape prediction errors were based on 
the actual prediction errors at the even numbers of strain- sensing stations. 

Analytical Surface Strains 

In the present paper, the surface strains, £■, (i = 0,l,2,3,....,n) , for inputs to the Displacement Transfer 
Functions were calculated analytically (not actually measured) by using the SPAR program. Namely, finite 
elements were used to simulate actual strain gages. The surface bending strains, (i = 0,l,2,3,....,n) , at the 
/-th strain-sensing stations (fig. 1) were generated by converting the nodal or element axial stresses into 
axial strains, (bending strains), through stress-strain relationship. 
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It must be mentioned that the surface strain, £. , can also be obtained from the span-wise nodal 
displacement differentials of the SPAR elements. However, it was found that in the beam-tip regions of the 
highly bent tapered beams (large deformations), the axial strains became negative because the output nodal 
displacements were the projected displacements along the original beam axis, and not along the deformed 
beam axis to reflect true span-wise displacements. Therefore, this displacement method was not used. 

Figure 5 shows the plots of SPAR-generated surface strains for the two loading cases, P - {300, 600} 
lb. Note that the strain level doubled when the applied load was increased from P - 300 lb to P = 600 lb 
because of linear elasticity used in the SPAR program. Note that the surface strains for the two loading 
cases increased practically linearly from the fixed ends and reached their respective peaks in the outboard 
region, and then decreased rapidly toward the beam tip (because of decreasing depth factor), and finally 
reached small nonzero values at the beam tip. Theoretically the beam-tip strains should be zero. Because 
of finite length, the beam-tip finite element (simulated strain gage), which has one end at the beam tip, but 
the other end is slightly off from the beam tip, gave small nonzero strains. By reducing the element sizes, 
the small nonzero beam-tip strains could be reduced. Keep in mind that if an actual strain gage is installed 
at the beam tip, the strain output could also be nonzero because of finite length of the strain gage. 

Criteria of Prediction Errors 

The SPAR-calculated surface strain data of figure 5 were then input to the Displacement Transfer 
Functions (constant-domain and variable-domain cases) to calculate the slopes and deflections. The shape 
prediction errors of the Displacement Transfer Functions were estimated by comparison with the SPAR- 
calculated slopes and deflections, which were used as reference yardsticks. The prediction errors are then 
defined as the percent of normalized differences (signs neglected) between the predicted slopes and 
deflections, {(^)pred >(T)pred 1 ’ corresponding SPAR-generated slopes and deflections, 

m SPAR. ’(Ji) spar} • Namely, 


Slope ( 6*, ) prediction error = 


Deflection ( j; ) prediction error = 


(^;)pred. (^i)sPAR 

X 100% 

(i = 

1,2,3,.. 

..,n) 

(23) 

(^i)sPAR 

( 3^1 )pi-ed. ~ ( 3 ^; )sPAR 

X 100% 

(i = 

1,2,3,.. 

..,n) 

(24) 

(^lOsPAR 


SELECTED DISPLACEMENT TRANSFER FUNCTIONS 

Because the cantilever tubular beam is highly tapered, the following four Displacement 
Transfer Functions for nonuniform beams were chosen in the prediction accuracy analysis. 

Based on piecewise-linear strain representations: 

1. Equations (5) — constant-domain Nonuniform Displacement Transfer Functions formulated for 
nonuniform beams based on piecewise-linear strain representations. 

2. Equations (13) — variable-domain Nonuniform Displacement Transfer Eunctions formulated for 
nonuniform beams based on piecewise-linear strain representations. 

Based on piecewise-nonlinear strain representations: 

3. Equations (7) — constant-domain Improved Displacement Transfer Eunctions formulated for 
nonuniform beams based on piecewise-nonlinear strain representations. 

4. Equations (17) — variable-domain Improved Displacement Transfer Eunctions formulated for 
nonuniform beams based on piecewise-nonlinear strain representations. 
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The above four Displacement Transfer Functions were used to calculate both slopes and deflections for 
comparison with the corresponding SPAR-calculated values. 

NUMERICAL RESULTS 

The SPAR-generated surface strain data of figure 5 were then input into equations (5a), (5b), (5c), 
(13a), (13b), and (13c) (piecewise-linear strain representations), and equations (7a), (7b), (7c), (17a), (17b), 
and (17c) (piecewise-nonlinear strain representations) for the calculations of both slopes and deflections. 
The slope and deflection prediction errors were then calculated respectively from the prediction error 
equations (23), and (24). 


Calculated Slope Data 

The slopes, 0 - , calculated from the slope equations (5a), and (13a) (based on piecewise-linear strain 
representations) and slope equations (7a), and (17a) (based on piecewise-nonlinear strain representations 
for the two loading cases, P - {300, 600} lb, are listed below. The percent slope prediction errors were 
calculated from the slope prediction error equation (23). 

For Piecewise-Linear Strain Representations: 

The slope prediction errors based on piecewise-linear strain representations are listed in tables 2 and 3 
for loading cases of P = (300, 600} lb. 
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Table 2. Comparisons of slopes, , calculated from SPAR and from constant and variable-domain slope 
equations [(5a), (13a)] (piecewise-linear strain representations); tapered cantilever tubular beam 
( I = 300 in., Cg = 4 in., = 1 in.) subjected to tip load of P = 3001b; n = 16 , A/ = 18.75 in., fine 
(A/), = 18.75 in., coarse (A/)- =37.50 in. 


i 

SPAR 
6i , deg 

Constant-domain [eq. (5a)] 

01 , deg Percent error 

Variable-domains [eq. (13a)] 
6j , deg Percent error 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1 

2.1815 

2.0895 

4.2174 

2.0895 

4.2174 

2 

4.3526 

4.3460 

0.1502 

4.3460 

0.1502 

3 

7.0977 

6.7873 

4.3742 

6.7873 

4.3742 

4 

9.4255 

9.4213 

0.0447 



5 

12.3764 

12.2534 

0.9943 

12.2542 

0.9876 

6 

15.2999 

15.2918 

0.0525 

15.2926 

0.0473 

7 

18.9492 

18.5395 

2.1623 

18.5403 

2.1582 

8 

21.9903 

21.9831 

0.0326 



9 

25.1656 

25.5988 

1.7214 

25.5937 

1.7013 

10 

29.3683 

29.3544 

0.0474 

29.3497 

0.0634 

11 

33.0596 

33.1950 

0.4097 

33.1907 

0.3966 

12 

37.0656 

37.0367 

0.0780 



13 

40.9426 

40.7621 

0.4409 

40.6960 

0.6024 

14 

44.2216 

44.1650 

0.1278 

44.1057 

0.2620 

15 

46.8926 

46.8405 

0.1112 

46.7865 

0.2262 

16 

48.2679 

48.0961 

0.3559 

48.0446 

0.4625* 


30-percent increase from constant-domain case. 
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Table 3. Comparisons of slopes, , calculated from SPAR and from constant and variable-domain 
non-uniform slope equations [(5a), (13a)] (piecewise-linear strain representations); long tapered cantilever 
tubular beam (Z = 300 in., Cq = 4 in., = lin.) subjected to tip load of P = 6001b, n = 16, Al = 18.75 in.; 
fine (A/); = 18.75 in., coarse (A/), =37.50 in. 


SPAR Constant-domain [eq. (5a)] Variable-domains [eq. (13a)] 


i 

Oi , deg 

Oi , deg 

Percent error 

Of , deg 

Percent error 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1 

4.0121 

4.1735 

4.0242 

4.1735 

4.0242 

2 

8.6553 

8.6426 

0.1463 

8.6426 

0.1463 

3 

13.5897 

13.3892 

1.4750 

13.3892 

1.4750 

4 

18.3669 

18.3591 

0.0424 



5 

23.2829 

23.4785 

0.8400 

23.4800 

0.8462 

6 

28.6842 

28.6712 

0.0453 

28.6726 

0.0407 

7 

34.0679 

33.8507 

0.6375 

33.8519 

0.6340 

8 

38.9263 

38.9162 

0.0260 



9 

43.6006 

43.7767 

0.4039 

43.7703 

0.3891 

10 

48.3787 

48.3626 

0.0332 

48.3571 

0.0446 

11 

52.7958 

52.6120 

0.3481 

52.6075 

0.3567 

12 

56.4976 

56.4700 

0.0490 



13 

59.7895 

59.8850 

0.1597 

59.8269 

0.0626 

14 

62.8071 

62.7611 

0.0732 

62.7127 

0.1502 

15 

65.0539 

64.8797 

0.2679 

64.8381 

0.3318 

16 

65.9637 

65.8349 

0.1954 

65.7962 

0.2540* 


30-percent increase from constant-domain case. 


Note from tables 2, and 3 (piecewise-linear strain representations) that, by removing three strain- 
sensing stations at i - 4, 8, 12 to change the constant-domain cases into variable-domain cases, the beam- 
tip slope (6*,) prediction errors increased from 0.3559 percent to 0.4625 percent (30-percent increase) for 
loading case of P = 300 lb (table 2), and from 0. 1954 percent to 0.2540 percent (also a 30-percent increase) 
for loading cases of P = 600 lb (table 3). However, the slope prediction errors are still in the negligible 
levels. Note also that by increasing the load from P = 300 lb to P = 600 lb, the beam-tip slope ( ) 
prediction errors decreased from 0.3559 percent (table 2) down to 0. 1954% (table 3) (45-percent reduction) 
for the constant-domain cases, and decreased from 0.4625 percent (table 2) down to 0.2540 percent (table 
3) (also a 45-percent reduction) for the variable-domain cases. 

For Piecewise-Nonlinear Strain Representations: 

The slope prediction errors based on piecewise-nonlinear strain representations are listed in tables 4 
and 5 for loading cases of P = (300, 600} lb. 
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Table 4. Comparisons of slopes, , calculated from SPAR and from constant and variable-domain slope 
equations [(7a), (17a)] (piecewise-nonlinear strain representations); tapered cantilever tubular beam 
( I = 300 in., Cq = 4 in., = 1 in.) subjected to tip load of T’ = 300lb; n = 16 , Al = 18.75 in., fine 
(A/), = 18.75 in., coarse (A/)- =37.50 in. 


SPAR Constant-domain [eq. (7a)] Variable-domains [eq. (17a)] 


i 

0i , deg 

0i , deg 

Percent error 

, deg 

Percent error 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1 

2.1815 

2.0884 

4.2683 

2.0884 

4.2683 

2 

4.3526 

4.3447 

0.1813 

4.3447 

0.1813 

3 

7.0977 

6.7862 

4.3890 

6.7862 

4.3890 

4 

9.4255 

9.4201 

0.0571 



5 

12.3764 

12.2516 

1.0089 

12.2655 

0.8964 

6 

15.2999 

15.2903 

0.0629 

15.3038 

0.0258 

7 

18.9492 

18.5391 

2.1644 

18.5522 

2.0953 

8 

21.9903 

21.9837 

0.0300 



9 

25.1656 

25.6010 

1.7301 

25.6218 

1.8130 

10 

29.3683 

29.3606 

0.0262 

29.3801 

0.0400 

11 

33.0596 

33.2072 

0.4467 

33.2252 

0.5010 

12 

37.0656 

37.0587 

0.0184 



13 

40.9426 

40.8057 

0.3343 

40.7581 

0.4506 

14 

44.2216 

44.2497 

0.0636 

44.2071 

0.0328 

15 

46.8926 

46.9762 

0.1784 

46.9376 

0.0959 

16 

48.2679 

48.2910 

0.0479 

48.2543 

0.0283* 


41 -percent reduction from constant-domain case. 
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Table 5. Comparisons of slopes, , calculated from SPAR and from constant and variable-domain slope 
equations [(7a), (17a)] (piecewise-nonlinear strain representations); tapered cantilever tubular beam 
( I = 300 in., Cq=4 in., =1 in.) subjected to tip load of P = 600 lb; n = 16 , Al = 18.75 in., fine 
(AZ)j = 18.75 in., coarse (AZ); = 37.50 in. 


i 

SPAR 
di , deg 

Constant-domain [eq. (7 a)] 

9j , deg Percent error 

Variable-domains [eq. (17a)] 
9i , deg Percent error 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1 

4.0121 

4.1713 

3.9690 

4.1713 

3.9690 

2 

8.6553 

8.6400 

0.1771 

8.6400 

0.1771 

3 

13.5897 

13.3872 

1.4898 

13.3872 

1.4898 

4 

18.3669 

18.3569 

0.0542 



5 

23.2829 

23.4753 

0.8263 

23.4999 

0.9317 

6 

28.6842 

28.6686 

0.0545 

28.6911 

0.0238 

7 

34.0679 

33.8501 

0.6393 

33.8702 

0.5803 

8 

38.9263 

38.9170 

0.0238 



9 

43.6006 

43.7796 

0.4104 

43.8063 

0.4716 

10 

48.3787 

48.3698 

0.0183 

48.3924 

0.0284 

11 

52.7958 

52.6249 

0.3237 

52.6438 

0.2879 

12 

56.4976 

56.4911 

0.0116 



13 

59.7895 

59.9232 

0.2237 

59.8815 

0.1538 

14 

62.8071 

62.8299 

0.0364 

62.7953 

0.0188 

15 

65.0539 

64.9841 

0.1073 

64.9544 

0.1530 

16 

65.9637 

65.9811 

0.0263 

65.9536 

0.0154* 


41 -percent deduction from constant-domain case. 


Note from tables 4, and 5 (piecewise-nonlinear strain representations) that, by removing three strain- 
sensing stations at i - 4, 8, 12 to change the constant-domain cases into variable-domain cases, the beam- 
tip slope (6(j) prediction errors decreased from 0.0479 percent down to 0.0283 percent (41 -percent 
reduction) for the loading case of P = 300 lb (table 4), and from 0.0263 percent down to 0.0154 percent 
(also a 41 -percent reduction) for loading case of P = 600 lb (table 5) [remember, 30-percent increase for 
stepwise linear strain representations (tables 2, and 3)]. 

Also note that by increasing the load from P = 300 lb to P = 600 lb, the beam-tip slope ( ) prediction 
errors decreased from 0.0479 percent (table 4) down to 0.0263 percent (table 5) (45-percent reduction) for 
the constant-domain cases, and decreased from 0.0283 percent (table 4) down to 0.0154 percent (table 5) 
(also a 45-percent reduction) for the variable-domain cases. 

It is important to mention that by changing the piecewise-linear strain representations (tables 2, and 3) 
into piecewise-nonlinear strain representations (tables 4, and 5), the beam-tip slope ( ) prediction errors 
for both loading cases, P = (300, 600} lb, were reduced by as large as 87 percent for the constant-domain 
cases, and by 94 percent for the variable-domain cases. Note from tables 4 and 5 that nonlinear strain 
representations reducing the number of domains (that is, increasing domain lengths) could somewhat 
improve the beam-tip slope (6*„) prediction accuracies. Such domain-length-dependent shape prediction 
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accuracies of the Displacement Transfer Functions were also observed in the experimental shape-prediction 
accuracy studies (ref. 13). 


Calculated Slope Curves 

Figures 6 and 7 respectively show the slope ( 6*, ) curves generated using the slope data of tables 2 and 
3 (piecewise-linear strain representations), and tables 4 and 5 (piecewise-nonlinear strain representations). 
Note from figures 6 and 7 that the slope-curves are slightly s-shaped, and the values of 6*, did not double 
when the load was increased from P = 300 lb to F = 600 lb. Namely, 6*, is not a linear function of P. 
Because of infinitesimal slope prediction errors, the predicted slope curves for all the constant and variable- 
domain cases pictorially fell on top of the corresponding SPAR slope curves, showing remarkable 
accuracies of the Displacement Transfer Functions formulated based on piecewise approaches. 

Calculated Deflection Data 

The deflections, y - , calculated from the deflection equations (5b), (5c), (13b), and (13c) (based on 
piecewise-linear strain representations) and deflection equations (7b), (7c), (17b), and (17c) (based on 
piecewise-nonlinear strain representations) for the two loading cases, P - {300, 600} lb are listed below. 
The percent deflection prediction errors were calculated from the deflection prediction error equation (24). 

For Piecewise-Linear Strain Representations: 

The deflections, y^, calculated from the deflection equations (5b), (5c) (13b), and (13c)] (based on 
piecewise-linear strain representations) for the loading cases, P - (300, 600} lb are listed respectively in 
tables 6, and 7. 
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Table 6. Comparisons of deflections, , calculated from SPAR and from constant and variable-domain 
deflection equations [(5b, 5c), (13b, 13c)] (piecewise-linear strain representations); tapered cantilever 
tubular beam (/ = 300 in., Cq =4in., =lin.) subjected to tip load of P = 3001b, « = 16 , A/ = 18.75 in., 
fine (AZ)i = 18.75 in., coarse (AZ)j = 37.50 in. 


i 

SPAR 
Ji , in. 

Constant domain [eqs. (5b,5c)] 
y, Percent error 

Variable domains [eqs. (13b, 13c)] 
y, Percent error 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1 

0.3223 

0.3377 

4.7845 

0.3377 

4.7845 

2 

1.4029 

1.3871 

1.1270 

1.3871 

1.1269 

3 

3.3148 

3.2096 

3.1741 

3.2096 

3.1741 

4 

5.8910 

5.8746 

0.2789 



5 

9.3147 

9.4591 

1.5503 

9.4594 

1.5533 

6 

14.0737 

14.0501 

0.1678 

14.0506 

0.1638 

7 

20.0231 

19.7479 

1.3747 

19.7487 

1.3705 

8 

26.6907 

26.6659 

0.0928 



9 

34.5995 

34.9302 

0.9558 

34.9294 

0.9533 

10 

44.7190 

44.6814 

0.0841 

44.6785 

0.0906 

11 

56.6063 

56.0742 

0.9399 

56.0693 

0.9486 

12 

69.3198 

69.2690 

0.0734 



13 

83.8382 

84.4150 

0.6880 

84.3710 

0.6355 

14 

101.7093 

101.6065 

0.1010 

101.5248 

0.1813 

15 

121.7164 

120.7512 

0.7930 

120.6318 

0.8910 

16 

141.6016 

141.3056 

0.2090 

141.1485 

0.3199* 


* 53-percent increase from constant-domain case. 
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Table 7. Comparisons of deflections, T, , calculated from SPAR and from constant and variable-domain 
deflection equations [(5b, 5c), (13b, 13c)] (piecewise-linear strain representations); tapered cantilever 
tubular beam ( Z = 300 in., CQ=4in., =1 in.) subjected to tip load of P = 600 lb, n = 16, A/ = 18.75in., 
fine (AZ)i = 18.75 in., coarse (AZ)j = 37.50 in. 


i 

SPAR 
T/ , in. 

Constant domain [eqs. (5b, 5c)] 
y, Percent error 

Variable domains [eqs. (13b, 13c)] 
y, Percent error 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1 

0.6446 

0.6755 

4.7845 

0.6755 

4.7845 

2 

2.8058 

2.7742 

1.1270 

2.7742 

1.1270 

3 

6.6296 

6.4192 

3.1741 

6.4192 

3.1741 

4 

11.7821 

11.7492 

0.2789 



5 

18.6294 

18.9182 

1.5503 

18.9187 

1.5533 

6 

28.1474 

28.1001 

0.1678 

28.1013 

0.1638 

7 

40.0463 

39.4957 

1.3747 

39.4974 

1.3705 

8 

53.3814 

53.3319 

0.0928 



9 

69.1991 

69.8605 

0.9558 

69.8587 

0.9533 

10 

89.4381 

89.3628 

0.0841 

89.3570 

0.0906 

11 

113.2126 

112.1485 

0.9399 

112.1386 

0.9486 

12 

138.6397 

138.5379 

0.0734 



13 

167.6765 

168.8301 

0.6880 

168.7421 

0.6355 

14 

203.4186 

203.2131 

0.1010 

203.0497 

0.1813 

15 

243.4327 

241.5024 

0.7930 

241.2637 

0.8910 

16 

283.2032 

282.6112 

0.2090 

282.2971 

0.3199* 


* 53-percent increase from constant-domain case. 


Note from tables 6, and 7 (piecewise-linear strain representations) that deflections ( y-) were doubled when 
tbe load was doubled from P = 300 lb to P = 600 lb because y- are linear functions of P. Note also that 
tbe _y. -prediction percent errors under tbe two loading cases, P = (300, 600} lb, (tables 6, and 7) are 

identical for both constant and variable-domain cases, and are independent of tbe loading magnitude of P. 

As shown in tables 6, and 7, by removing three strain-sensing stations at i - 4, 8, 12, to change the 
constant-domain cases to variable-domain cases, the beam-tip deflection ( _y„ ) prediction errors increased 

from 0.2090 percent (constant-domain cases) to 0.3199 percent (variable-domain cases) (53-percent 
increase) for both loading cases, P = (300, 600} lb. 

For Piecewise-Nonlinear Strain Representations: 

The deflections, y, , calculated from the deflection equations (7b), (7c) (17b), and (17c) (based on 
piecewise-nonlinear strain representations) for the loading cases, P = (300, 600} lb are listed respectively 
in tables 8, and 9. 
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Table 8. Comparisons of deflections, y. , calculated from SPAR and from constant and variable-domain 
deflection equations [(7b, 7c), (17b, 17c)] (piecewise-nonlinear strain representations); tapered cantilever 
tubular beam ( I = 300 in., Cq=4 in., = 1 in.) subjected to tip load of 7* = 300 lb, n - 16, fine 
(A/); = 18.75 in., coarse (AZ), =37.50 in. 


i 

SPAR 
y . , in. 

Constant domain [eqs. (7b,7c)] 
y - , in. Percent error 

Variable domains [eqs. (17b, 17c)] 
y ^ , in. Percent error 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1 

0.3223 

0.3376 

4.7284 

0.3376 

4.7284 

2 

1.4029 

1.3865 

1.1688 

1.3865 

1.1688 

3 

3.3148 

3.2086 

3.2038 

3.2086 

3.2038 

4 

5.8910 

5.8732 

0.3019 



5 

9.3147 

9.4572 

1.5303 

9.4619 

1.5800 

6 

14.0737 

14.0476 

0.1853 

14.0570 

0.1184 

7 

20.0231 

19.7451 

1.3887 

19.7592 

1.3179 

8 

26.6907 

26.6632 

0.1032 



9 

34.5995 

34.9280 

0.9494 

34.9554 

1.0286 

10 

44.7190 

44.6810 

0.0851 

44.7168 

0.0051 

11 

56.6063 

56.0780 

0.9333 

56.1222 

0.8552 

12 

69.3198 

69.2812 

0.0557 



13 

83.8382 

84.4453 

0.7241 

84.4750 

0.7595 

14 

101.7093 

101.6760 

0.0327 

101.6786 

0.0302 

15 

121.7164 

120.8950 

0.6749 

120.8703 

0.6951 

16 

141.6016 

141.5683 

0.0235 

141.5165 

0.0601* 


156-percent increase from constant-domain case. 
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Table 9. Comparisons of deflections, y. , calculated from SPAR and from constant and variable-domain 
deflection equations [(7b, 7c), (17b, 17c)] (piecewise-nonlinear strain representations); tapered cantilever 
tubular beam ( 1 = 300 in., Cq= 4 in., c„ = 1 in.) subjected to tip load of T’ = 600lb, u = 16, fine 
(A/), = 18.75 in., coarse (A/)- =37.50 in. 


i 

SPAR 
y, , in. 

Constant domain [eqs. (7b,7c)] 
y - , in. Percent error 

Variable domains [eqs. (17b, 17c)] 
y - , in. Percent error 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1 

0.6446 

0.6751 

4.7282 

0.6751 

4.7282 

2 

2.8058 

2.7730 

1.1688 

2.7730 

1.1688 

3 

6.6296 

6.4172 

3.2038 

6.4172 

3.2038 

4 

11.7821 

11.7465 

0.3019 



5 

18.6294 

18.9145 

1.5303 

18.9237 

1.5800 

6 

28.1474 

28.0952 

0.1853 

28.1140 

0.1184 

7 

40.0463 

39.4901 

1.3887 

39.5185 

1.3179 

8 

53.3814 

53.3264 

0.1032 



9 

69.1991 

69.8560 

0.9494 

69.9109 

1.0286 

10 

89.4381 

89.3619 

0.0851 

89.4335 

0.0051 

11 

113.2126 

112.1560 

0.9333 

112.2443 

0.8552 

12 

138.6397 

138.5624 

0.0557 



13 

167.6765 

168.8907 

0.7241 

168.9500 

0.7595 

14 

203.4186 

203.3521 

0.0327 

203.3571 

0.0302 

15 

243.4327 

241.7899 

0.6749 

241.7406 

0.6951 

16 

283.2032 

283.1367 

0.0235 

283.0330 

0.0601* 


* 156-percent increase from constant-domain case. 


Note from tables 8 and 9 (piecewise-nonlinear strain representations) that deflections ( ) were doubled 

when the load was doubled from P - 300 lb to P = 600 lb because are linear functions of P. It is 
important to note that the _y^. -prediction errors are identical under the two loading cases, P - (300, 600} lb, 

(tables 8 and 9) for both constant and variable-domain cases and are, therefore, invariant to the loading 
magnitude of P. 

By removing three strain- sensing stations at i - 4, 8, 12, to convert the constant-domain cases to 
variable-domain cases, the beam-tip deflection ( y„ ) prediction errors increased from 0.0235 percent 
(constant-domain cases) to 0.0601 percent (156-percent increase) for both loading cases, 
P- (300, 600} lb. [Remember, there is only a 53-percent increase for piecewise-linear strain 
representations (tables 6, and 7)]. 


Calculated Deflection Curves 

Figures 8 and 9 respectively show the deflection ( y . ) curves generated using the deflection data of 
tables 6, and 7 (piecewise-linear strain representations) and tables 8, and 9 (piecewise-nonlinear strain 
representations). Note from figures 8 and 9 that the y, -curves are shallow bow-shaped, and the values of 
y. were doubled when the load was increased from P - 300 lb to P = 600 lb. Namely, y. is a linear 
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function of P. Because of infinitesimal deflection prediction errors, the predicted deflection curves for all 
the constant and variahle-domain cases practically fell on top of the corresponding SPAR deflection curves, 
showing amazing high accuracies of the Displacement Transfer Functions formulated based on piecewise 
formulations. 


SHAPE PREDICTION ERROR ANALYSIS 

This section compares the shape prediction errors of constant and variahle-domain Displacement 
Transfer Functions formulated using linear or nonlinear strain representations. The slope and deflection 
prediction errors were calculated respectively from the shape prediction error equations (23) and (24). 

Slope Prediction Error Curves 

Figures 10 and 11 respectively show the slope (6*, ) prediction error curves for loading cases of 
P - {300, 600} Ih based on the slope data listed in tables 2 and 3 (piecewise-linear strain representations). 
The beam-tip slope {0^) prediction errors (in percent) are also indicated in figures 10 and 11. Those 
0^ -error curves were plotted by using only the actual prediction error data at the even strain-sensing stations 
which are either coincidental with SPAR nodes or at the midpoints between two adjacent SPAR axial nodes. 
The larger prediction errors at the odd strain-sensing stations (tables 2 and 3) cannot represent the true 
prediction errors because the odd strain-sensing stations are not precisely at the SPAR nodes and, therefore, 
were not plotted in figures 10 and 11. 

The constant-domain cases and variable-domain cases are shown respectively with solid curves (with 
solid circular symbols), and dashed curves (with open square symbols). The larger errors near the inboard 
regions are due to small numbers divided by other small numbers. The variable-domain cases have slightly 
larger errors in the outboard regions (figs. 10 and 11). 

Note from figures 10 and 1 1 that, doubling the load from P = 300 lb to P = 600 lb, the 6*, -prediction 
errors in the outboard regions came down considerably (45 -percent error reduction at the beam tip for both 
constant-domain and variable-domain cases). 

Figures 12 and 13 respectively show the predicted slope ( 6*, ) error curves for loading cases of 
P- (300, 600} lb based on the slope data listed in tables 4 and 5 (piecewise-nonlinear strain 
representations). The constant-domain cases and variable-domain cases are shown respectively with solid 
curves (with solid circular symbols), and dashed curves (with open square symbols). The amounts of the 
beam- tip slope ( ) prediction errors are also shown in figures 12 and 13. 

Note from figures 12 and 13 that, doubling the load from P = 300 lb to P = 600 lb, the 6*, -prediction 
errors in the outboard regions came down considerably, and at the beam tip, the prediction errors also 
decreased by 45 percent for both constant-domain and variable-domain cases. 

It is very important to mention that by changing the piecewise-linear strain representations (figs. 10 and 
11) into piecewise-nonlinear strain representations (figs. 12 and 13), the beam-tip slope (6},) prediction 
errors for both loading cases, P = (300, 600} lb could be reduced remarkably. Namely, reducing from 
(0.3559 percent, 0.1954 percent} (figs. 10 and 1 1) down respectively to (0.0479 percent, 0.0263 percent} 
(87-percent error reduction) (figs. 12 and 13) for the constant-domain cases, and reducing from 
(0.4625 percent, 0.2540 percent} (figs. 10 and 1 1) down respectively to (0.0283 percent, 0.0154 percent} 
(94-percent error reduction) (figs. 12 and 13) for the variable-domain cases. 
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Deflection Prediction Error Curves 

Figures 14 and 15 respectively show the predicted deflection ( y - ) error curves based on the deflection 
error data at even strain-sensing stations listed in tables 6 and 7 (piecewise-linear strain representations) 
and tables 8 and 9 (piecewise-nonlinear strain representations). The predicted deflection ( y, ) error curves 
are identical for both loading cases, P - {300, 600} lb. The constant-domain cases and variable-domain 
cases are shown respectively with solid curves (with solid circular symbols), and dashed curves (with open 
square symbols). Again, the large errors (figs. 14 and 15) near the fixed end are due to a small number 
divided by another small number. 

As shown in figures 14 and 15, changing from piecewise-linear strain representations into 
piecewise-nonlinear strain representations, the y, -prediction errors could be reduced considerably in the 
outboard region. For example, the beam-tip deflections ( y„ ) prediction errors for both loading cases, 
P - (300, 600} lb could be reduced from 0.2090 percent down to 0.0235 percent (89 percent error 
reduction) for the constant-domain cases, and from 0.3199 percent down to 0.0601 percent(81 percent error 
deduction) for the variable-domain cases. 

EXPERIMENTAL VALIDATIONS OF SHAPE-PREDICTION 

ACCURACIES 

In the process of writing the present report, large-scale ground loads tests of the Global Observer wing 
( 175-ft. wingspan) were carried out at NASA Dry den (now Armstrong) Flight Research Center Flight Loads 
Laboratory (Edwards, California) (ref. 14). The experimentally measured strains obtained from the ground 
tests were then used to validate the Displacement Transfer Functions accuracies. The Global Observer wing 
is slightly tapered and consisted of several sections joined together. Therefore, the strain-sensing stations 
had to be properly distributed in order to avoid positioning them at the connection junctures, and thus 
creating the variable-domain case. The surface strains were measured by using four fiber-optic 
strain- sensing lines, with two lines on the wing upper surface, and another two lines on the wing lower 
surface. The whole wing was loaded from 0 percent up to 100 percent Design Limit Load (DLL) at which 
the wing-tip slope exceeded 20 degrees, and wing-tip deflection reached nearly 14 ft. The variable-domain 
Displacement Transfer Function (14c), developed for the slightly nonuniform beam, was used for the wing 
shape calculations. Using the measured surface strain data, the wing deflections were calculated from 
equation (14c), and were compared with the photogrammetry data, which were used as yardsticks to 
estimate the shape prediction errors. 

Figures 16, and 17 (duplications of figures 11 and 12 of ref. 14) respectively show the plots of the 
predicted and measured wing deflections along the strain-sensing lines Number 3 and Number 4. Note that 
the strain data points are not uniformly distributed because of variable-domains. Amazingly, the measured 
data points practically fell on top of the corresponding predicted deflection curves for the whole range of 
loading levels (0-1 00-percent DLL). The wing-tip deflection prediction errors based on Number 3 and 
Number 4 strain-sensing lines are respectively in the small ranges of (0.09-1.72 percent) and 
(0.25-1.45 percent) at 100-percent DLL, depending on the strain-sensing domain distances used (ref. 14). 

In this ground test, it was found that by increasing the domain length from (A/); = 1 in. to 
(A/); = 12 in., the averaged wing-tip deflection prediction errors at 100-percent DLL could be reduced 
from 1.92 percent down to 0.27 percent (a remarkable 86-percent error reduction). This particular case 
certainly violates the principle that reducing the strain-sensing domains (that is, increasing the number of 
strain- sensing stations) can improve the shape prediction accuracies. Because of noisy strain data, 
increasing the number of strain-sensing stations caused more errors to be accumulated (fig. 10 of ref. 14). 
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This large-scale ground test results thus reinforced the previous finite-element analytical validation results 
(refs. 1-9), and give confidence in the shape prediction accuracies of the Variahle-Domains Displacement 
Transfer Functions. 

FACTORS AFFECTING SHAPE PREDICTION ACCURACIES 

The surprising high accuracies of analytical deformed shape predictions could he attributed to the 
piecewise (or recursive) formulations of the Displacement Transfer Functions. The major factors affecting 
the shape prediction accuracies of the Displacement Transfer Functions can he listed as: 

1. Domain lengths (that is, intervals of strain-sensing stations along each strain-sensing line), 

2. Accuracies of piecewise-linear representations of depth-factor distributions, 

3. Accuracies of piecewise-linear or piecewise-nonlinear representations of strain distributions, 

4. Accuracies of the input surface strains, 

5. Accuracies of depth factors determined from surface strains. 

CONCLUDING REMARKS 

The variable-domain Displacement Transfer Functions were formulated through piecewise integrations 
of the shifted curvature-strain differential equation. Using a long tapered cantilever tubular beam, the shape 
prediction accuracies of the variable-domain Displacement Transfer Functions were compared with the 
shape prediction accuracies of the constant-domain Displacement Transfer Functions formulated earlier. 

Finite-element analysis was used to analytically generate surface strains. Also, the finite-element 
calculated slopes and deflections were used as reference yardsticks in the shape prediction error estimations. 
Some key results are listed below: 

1. The slope and deflection equations in recursive forms developed for the constant-domain cases could 
be converted to those for the variable-domain cases by simply replacing the constant-domain lengths with 
variable-domain lengths. 

2. The variable-domain Displacement Transfer Functions are highly accurate like the constant-domain 
Displacement Transfer Functions in shape predictions. 

3. The shape prediction errors (slopes and deflections) are in the negligible ranges, with beam-tip slope 
( ) prediction errors in the ranges of {(0.0263-0.3559 percent) and (0.0154-0.4625 percent)} respectively 
for the constant-domain and variable-domain cases; and the beam-tip deflection ( y„ ) prediction errors in 
the ranges of {(0.0253-0.2090 percent) and (0.0601-0.3199 percent)} respectively for the constant-domain 
and variable-domain cases. 

4. Based on either piecewise-linear or piecewise-nonlinear strain representations, the deflection 
prediction errors of the constant-domain and variable-domain Displacement Transfer Functions are 
invariant to the loading magnitudes. However, the slope prediction errors decreased with increasing loading 
magnitude. 

5. For piecewise-linear strain representations, changing from constant-domain cases into variable- 
domain cases, the beam-tip slope and deflection {0„, y,J prediction errors increased by {30 percent and 
53 percent} respectively, regardless of loading magnitude. 

6. For piecewise-nonlinear strain representations, changing from constant-domain cases into variable- 
domain cases, regardless of loading magnitude, the beam-tip slope (6*„) prediction errors decreased by 
41 percent, however, the beam-tip deflection (y„) prediction errors increased by 156 percent. 

7. By changing from the piecewise-linear strain representations to piecewise-nonlinear strain 
representations, the beam-tip slope ( ) prediction errors could be reduced by 87 percent for the constant- 
domain cases, and reduced by 94 percent for the variable-domain cases. 
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8. By changing from the piecewise-linear strain representations to piecewise-nonlinear strain 
representations, the heam-tip deflection ( y,, ) prediction errors could he reduced hy 89 percent for the 
constant-domain cases, and reduced hy 81 percent for the variahle-domain cases. 
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FIGURES 



Figure 1 . Embedded beam (cross section of structure, including a tubular beam) along strain-sensing line 
with evenly distributed strain-sensing stations. 



Piecewise linear strain representation: 

E(X) = E/_-| + (E/- X £ Xj) 

Piecew ise nonlinear strain representation: 


E(x) = E (A/)/+i [2(A/)/ + (A/)/+i]E/_-| — [(A/)/ + (A/),y-|] E/ + (A/)y E /^.1 ^ % 

(A/),(A/)y^i [(A/), + (A/),v,i] ' 


I (A/)/4.iE/_i — [(A/)/ + (A/)/+i]e/ + (A/)/E/+t 

(A/),(A/),-^i [(A/), + (A/),vi] 


(X,_i<X<X/) 


Figure 2. Piecewise linear and piecewise nonlinear representations of variable domain surface strain 
distribution. 
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Figure 3. Segment of deformed embedded beam used for geometrically relate lower and upper surface 
strains, {s., €■}, to lower and upper depth factors, {c., C-} . 


Size of SPAR model 



Figure 4. Large deformations of tapered cantilever tubular beam (SPAR finite-element model) subjected to 
upward tip load of P - 300 lb (or P - 600 lb); / = 300 in., Cq = 4 in., = 1 in. 
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Figure 5. SPAR-generated surface strain curves for tapered cantilever tubular beam subjected to to upward 
tip load of P = 300 lb (or P - 600 lb); I - 300 in., Cg = 4 in., c^- 1 in. 



Strain-sensing station, / 

Figure 6. Slopes, 0j, of tapered cantilever tubular beam calculated from constant domain and variable 
domain slope equations [(5a), and (13a)] (piecewise linear strain representations); I - 300 in., Cq - 4 in., 
= 1 in.;P= (300, 600} lb. 
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Figure 7. Slopes, 9i, of tapered cantilever tubular beam calculated from constant domain and variable 
domain slope equations [(7a), and (17a)] (piecewise nonlinear strain representations); / = 300 in., Cq -4 in., 
c^- I in.; P - (300, 600} lb. 



Strain-sensing station, / 

Figure 8. Deflections, y^, of tapered cantilever tubular beam calculated from constant domain and variable 
domain deflection equations [(5b, (5c), (13b), and (13c)] (piecewise linear strain representations); 
/ = 300 in., Cq-4 in., = 1 in.; P - (300, 600} lb. 
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Figure 9. Deflections, y^, of tapered cantilever tubular beam calculated from constant domain and variable 
domain deflection equations [(7b), (7c), (17b), and (17c)] (piecewise nonlinear strain representations); 
/ = 300 in., Cq-A in., = 1 in.; P - (300, 600} lb. 



Strain-sensing station, / 


Figure 10. Comparison of slope (0j) prediction errors of constant domain and variable domain slope 
equations [(5a), and (13a)] (piecewise linear strain representations); tapered cantilever tubular beam; 
/ - 300 in, Cq -4 in., c^- I in.; P - 300 lb. 
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Figure 11. Comparison of slope (0;) prediction errors of constant domain and variable domain slope 
equations [(5a), and (13a)] (piecewise linear strain representations); tapered cantilever tubular beam; 
/ = 300 in, Cq-A in., c^-\ in.; P - 600 lb. 



Strain-sensing station, / 


Figure 12. Comparison of slope (0;) prediction errors of constant domain and variable domain slope 
equations [(7a), and (17a)] (piecewise nonlinear strain representations); tapered cantilever tubular beam; 
/ - 300 in, Cq -4 in., = 1 in.; P - 300 lb. 
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Figure 13. Comparison of slope (0;) prediction errors of constant domain and variable domain slope 
equations [(7a), and (17a)] (piecewise nonlinear strain representations); tapered cantilever tubular beam; 
/ = 300 in, Cq-A in., c^-\ in.; P - 600 lb. 



Strain-sensing station, / 


Figure 14. Comparisons of deflection (y^) prediction errors of constant domain and variable domain 
deflection equations [(5b), (5c), (13b), and (13c))] (piecewise linear strain representations); tapered 
cantilever tubular beam; I - 300 in, Cq-4 in., Cy^-l in.; P - (300, 600} lb. 
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0.0601 % 
0 . 0235 % 


Figure 15. Comparisons of deflection (y;) prediction errors of constant domain and variable domain 
deflection equations [(7b), (7c)), (17b), and (17c)] (piecewise linear strain representations); tapered 
cantilever tubular beam; I - 300 in, Cq-4 in., c^- I in.; P - (300, 600} lb. 
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Figure 16. Comparison of predicted and measured wing deflections along No. 3 strain-sensing line for 
different loading levels (duplication of figure 11, ref. 14). 
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Figure 17. Comparison of predicted and measured wing deflections along No. 4 strain-sensing line for 
different loading levels (duplication of figure 12, ref. 14). 
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APPENDIX A 

DERIVATION OF STRAIN FUNCTION FOR VARIABLE DOMAIN 
PIECEWISE-NONLINEAR STRAIN REPRESENTATIONS 


Let ^{= X - Xj_i ) be the loeal axial eoordinate measured from the strain-sensing station, x,_j (fig. 2), 
and let the distribution of strain, £(<^) , within the two adjaeent domains, x,_[ <x<x- and x-<x< x,.^j, be 
presented by the following nonlinear function (quadratic equation) (fig. 2): 

f(^) = fl + K+ef ; ^ = x-x,._i (Ai) 

where a, b, and e are coefficients, which are to be determined in terms of three strains, £•,^ 1 } 

respectively at three strain- sensing stations, {x,_i , x, , x,^i } . 

At = 0(i.e., at strain-sensing station, x,_j ), we have £'(0) = £•,_) , and equation (Al) gives: 

a ^ £•,._! (A2) 

At = (Al)- (that is, at strain-sensing station, x, ), we have £'[(A/)J = s - , then in view of equation (A2), 
equation (Al) can be written as equation (A3): 

£i = + b(Al). + e(Al)f (A3) 


which can be rewritten as equation (A4): 


b + e(Al)j 



(A4) 


At ^ = (A/). -I- (A/), (that is, at strain- sensing station, x,^j ), we have £’[(A/), -I- (A/);_,_J = £';^j . Then in 
view of equation (A2), equation (Al) can be written as equation (A5): 

^,■,1 = + bmi + (Al),,, ] + eV(Al\ + (Al),„ f (A5) 

which can be rewritten as equation (A6): 


b + e[(Al),+(Al),„]^ 


(A/),+(A/),,, 


(A6) 


Subtracting equation (A4) from (A6), one can eliminate b and obtain the expression for e as equation (A7): 


e(Al) =—^ £j__£hL 

^ (A/),+(A/)„, (Al), 

(A/), (£,,, -£,_,)-[(AZ), +(A/),„](g, -£,_,) 
(A/),[(A/),+(A/),,J 
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(A/),g,^, - - [(A/), + + [(A/), + (A/),,Jg,_, 

(A/),[(A/), + (A/),,J 


(A7) 


From equation (A7), e can be written as: 


(A/),,,g,_,-[(A/), + (A/),„]fi + (A/),g,,, 
(A/),(A/),,,[(A/), + (A/),,J 


(A8) 


Substitution of equation (A8) into equation (A4) yields equation (A9): 


g,. -g,._i (A/),._,ig,._i - [(A/),. +(A/),.^Jg,. +(A/),.g,._,i 

(A/), (A/),(A/),,i[(A/),+(A/).,J 


(A/),,i[(A/),+(A/),,,](g, 


1 


g,_i)-(A/),(A/),,ig,,i +(A/),[(A/), +(A/),,i]g, -(A/)fg,,i 
(A/),(AZ),,i[(A/),+(A/),,J 

■-(AZ),,i[(AZ), +(AZ),,Jg;_i +(AZ),,i[(AZ), +(AZ),,Jg," 


(AZ),(AZ),,i[(AZ),+(AZ),,J 


(AZ),(AZ),,ig,_i +(AZ),[(AZ),. +(AZ),,Jg, -(AZ)fg, 


(A9) 


After grouping terms, one obtains the expression for h as equation (AlO): 


(AZ),,i[2(AZ), +(AZ),,Jg,_i -[(AZ), +(AZ),,J^g, +(AZ)fg,, 
(AZ),(AZ),,,[(AZ).+(AZ),,,] 


(AlO) 


In view of equations (A2), (A8), and (AlO), equation (Al) can be written as equation (All): 


£{x) = g,_i 


(AZ)„,[2(AZ), + (AZ)„,]g,_. - [(AZ), + (AZ),,.]^g, + (AZ)f g„. 

(AZ),(AZ),,,[(AZ), + (AZ),,J 
^ (AZ)„.g,_,-[(AZ),+(AZ),,Jg,+(AZ),g,,^ 

(AZ),(AZ),,J(AZ), + (AZ),,J 


(All) 


which is equation (1 1) in the text. 


For the constant-domain case, (AZ),_i = (AZ), = - AZ , equation (All) can be simplified as 

equation (A 12): 


£{x) 


3g 


(-1 


'i-l 


- 4g,. -r g, 
2AZ 


i+l 


{x- 


''/-l 


^ g,-i - + g,+i 2 


(A12) 


Equation (A12) is identical to equation (4a) in the text. 
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APPENDIX B 

DERIVATION OE VARIABLE-DOMAIN DISPLACEMENT TRANSEER 
EUNCTIONS EOR NONUNIEORM EMBEDDED BEAMS 

(Based on Piecewise-Linear Strain Representations) 

The slope equation (13a) and the deflection equation (13h) for nonuniform beams, based on piecewise- 
linear representations of both depth factors and surface strains, are duplicated below as equations (Bl) and 
(B2). 

Slope equation: 


tan6»,. = (A/),. 




■Ci (c,.-i-c,.) 


+ tan 


(Bl) 


(i = 1,2,3,... .,n) 


Deflection equation: 


J,-(A/)? 


2(c,._i-c,.) (c,_i-c.) 


c,Tog,^ + (c,._i-c,.) 


> + y,._i+(A/),.tan6';_i (B2) 

(i = i,2,3,....,n) 


Writing out equation (B2) for different indices, /(= 1,2,3,...., n) , and considering the indicial 
relationships expressed in equations (Bl) and (B2), one obtains equation (B3): 

For i = 1: 




^(^0 ^ i ) (^0 ^ i ) 


qlog, — + (Co-q) 

Cn 


' + }'o + (AOi tan6»o 


(B3) 
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For / = 2 : 




= ml 


+ (A/)f 


Cl C2 

C1C2 C2C1 

2 (ci -C2) 

(Ci — C2) 

Cl -C2 

C1C2 -C2C1 

2 (ci -C2) 

(Ci — C2 ) 

Cq “Cl 

CflCi — CiCq 


Cl log, — + (Ci -C2) 


^2 log, — + (Ci -C 2 ) 


■ + + (AZ)2 tan^j 


l2(c, 


0 Oi) (Cq C[) 

lan^i 


Cl log, — + (Co -Cl) 


+ (A/)2 (A/), 


^0 ^1 , ^ 0^1 ^iCq c^ 


Cq Cj (Cq Cl) 

After grouping terms, equation (B4) becomes: 


+ tan Or. 


>^0 +(A/)i tan6»o 


(B4) 


J2 =(A0^ 


+ (AZ)^ 


^2 ^iCl C 2 C 1 


2(Ci-C2) (C 1 -C 2 )' 


Sq ^f\C] S\C^ 


0^1 ^1^0 


2(co-Ci) (Co-Ci)" 


C2 log, — + (Ci -C2) 


Cllog,— + (Co - Cl) 


+ (A/)2 (AZ), 


^0 ^1 I ^qCi CiCq Cl 

Cq — Cl (Cq — Cl) Cq 


■ + 3^0 + [(AZ)i + (AZ)2]tan^o 


(B5) 
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For / = 3 : 


3^3 


^2 ^3 _ ^2^3 ^3^2 

2(C2-C3) (Cj-Cj)^ 


C 3 l 0 g, — + (C 2 -C 3 ) 


+ J 2 + (^^03 tan^2 


2 I ^2 ^3 ^2^3 ^3^2 


= (^03 1 3 

^M 2 (C 2 -C 3 ) (C 2 -C 3 )' 


C3l0g. — + (^ 2 -^ 3 ) 


C 3 

C 2 

V2 


+s 


{Mf 


S\ S 2 


^ 1^2 ^ 2^1 


^'2(Ci-C2) {c^-c^f 


C2lOg, — + (Ci-C2) 


+(A/)f 


^0 ^1 


£qC^ S^Cq 


2(Co-Ci) (Cq-Cj) 

+(A/)2(A/)j 


Cilog, — + (C(,-Ci) 


^0 ^1 , ^0^1 ^1^0 £l_ 


■ + ■ 


(Cq C\ ) 


+ }^o+[(A/)j+(A/)2]tan6»o 


tan ^2 


+ (A/)3 < (A/)2 


^1 ^2 _j_ ^1^2 ^2^1 

^1 ~ ^2 (^1 ~ ^2 )” 


log,— 

c, 


+(A/y 


^0 _|_ £c£i__£iSl}Qg AL 


(Cq Cj ) 


ii 

Cn 


+ tan 


(B6) 


After grouping terms, equation (B6) becomes equation (B7): 


J 3 =(A/)2< 

^2 ^3 

^ 2*'3 ^ 3^2 

C 3 l 0 g, ^ +(C2 C3) 

L C2 

2(C2 -C3) 

(^2 ^^ 3 ^ 

+ (A/)^ 

^1 -^2 

^1^2 “ ^2^1 

C2l0g,— + (Ci -C2) 

2(ci -C2) 

(Ci -C2)^ 

+ (A/)f< 


^0^1 ~ ^1^0 

Cl log,-^^ + (C(i-Ci) 

Co J 

0^ 

1 

0 

(Cq -Ci)^ 




^1 ^2 , ^ 1^2 ^ 2^1 


- + 


, 32-log^ — 

^1 ^2 (^1 ^ 2 ) 


+ [(A/)3+(A/)2]j(A/)i 

^0 ^1 (^0 ^ 1 ) 

+ Ao + [(A/)i + (A/)2 + (A/)3]tan0Q 


^0 ^1 gpCi ^1^0 


log,^ 


(B7) 
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For i = 4: 


J 4 

^ 3-^4 

^304 - ^ 4 C 3 

C 4 l 0 g,— + (C 3 -C 4 ) 

L ^3 jJ 

2 (C 3 -C 4 ) 

(C 3 -C 4 )^ 

= (AZ)^ 

^3 - ^4 

- .C 4 C 3 

C 4 l 0 g,— + (C 3 -C 4 ) 1 

C 3 JJ 

2 (C 3 -C 4 ) 



■ + 3^3 + (AZ)4 tan^3 


(A/); 


— S’x 


^2^3 ^3^2 


[2(C2-C3) {C2-C^f 


Cslog,— +(C 2 -C 3 ) 


+ (AZ)^ 


+ (AZ)f 


6\ S') 


S^C^-S'^C^ 


2 (Ci-C 2 ) (C1-C2) 

6*0 “^1 


C2l0g, — + (Ci -C 2 ) 


^0<-i 


+ (AZ)3 (AZ) 


^(^0 ^ 1 ) (^0 ^ 1 ) 

s^-s^ fjC2-^2^1 


Cllog, — + (Co - Cl) 


- + 


C1-C2 (C1-C2) 


log,— 




0^1 ^1^0 ^ 


+ [(AZ)3+(AZ)2]j(AZ)i 

^0 ^1 (Cq Cj) 

.+ 3^0 + [(AZ)i + (AZ)2 + (AZ)3]tan^Q 


tan 6 -i 


+ (AZ)4 


(AZ)3 
+ (AZ) 


^2 ^3 ^2^3 ^3^2 

^2 ~ O3 (C2 — C3) 


Sc\ S\ ^0^1 ^1^0 


log,^ 


+ (AZ), 


-6*^ S)C') - S')C 


- + - 


Cl c.^ 


1^2 ^2^1 

(Ci — C 2 ) 


log,— 


- + - 


Co -Cl (Cq-Ci) 


log,^ 


+ tan On 


(B8) 


47 


After grouping terms, equation (B8) becomes equation (B9): 


^4 =(A/)^ 


+ (AZ) 




+ {M)l 


^3 ^4 

^ 3^4 - ^ 4 C 3 

2(C3 -C4) 


£2 — £■3 

^ 2^3 “ ^ 3^2 

2(C2 -C3) 

(^2 ~ ^3) 

^1 “ ^2 

^1^2 ~ ^2^1 

rq 

1 

(Ci ~ C2) 

1 

(<r 

1 

0 

(0 

£qC^ — f [Cq 

1 

0 

(Cq — C[) 


Crlog, — + (C3 -C4) 


Cslog, — + (C2 -C3) 


C2lOg, — + (Ci -C2) 


Cilog, — + (Co-Ci) 


+ (AZ)4 (A/)3 




- + - 


2^3 *^3^2 

(^2 ~ ^3) 


log.^ 


■[(AZ)4+(AZ)3] (AZ), 


s, - s,c^ - s^c 


Cl C-) 


- + - 


1^2 ^2^1 
(Ci — C2) 


log.^ 


+ - 


0^1 "^1^0 


+ [(AZ)4 + (AZ)3 + (AZ)2 K (AZ)j 

[ ^ Cq “ Cj (Cq “ C[ ) 

+ Jo [(AZ)[ + (AZ)2 + (AZ)3 + (AZ)4]tan0g 


Cq Cj ^ 1^1 


log.^ 


(B9) 


Based on the indicial progression patterns in equations (B3), (B5), (B7), and (B9), the generalized form 
of the deflection equation can be written in dual summation forms (with different summation limits) as: 


^, = X (A/)’ 


j=l 








2(c,_ 


-c.) (c 


j-i 


;^uiog.^+(c,.-cp 


”+yo 


Contributions from deflection terms 


+X I + •■••+(AZ),](AZ) . 


e, ,c 




(c^-.-cp C._, 


■ + 


I (AO, 


J=l 


tan0„ 


Contributions trom slope terms 


(BIO) 
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Equation (BIO) can be rewritten in a more compact form as: 







2(c,,-cJ (c.._,-c) 


Cyl0g.-^ + (S-l-S), 


-t-1 


^ + 3^0 


Contributions from deflection terms 


+!• 

j=i 


i (A/), 


k=J+l 


(AO, 




S-i-c, 


(Cm-cO 




■ + 


I (A/), 


J=1 


tan0„ 


(Bll) 


Contributions from slope terms 


(i = 1,2,3,...., «) 


Equation (B1 1) is identical to equation (13c) in the text, and is called the Variable-Domain Nonuniform 
Displacement Transfer Eunction in dual summation form for nonuniform embedded beams based on 
piece wise-linear strain representations. 

Eor constant-domain lengths, (A/)j = (A /)2 = (A /)3 = = Al , equation (Bll) degenerates into: 






£j-fj ^fj-i 


^ 2(c._i-c.) (Cj_i~Cj) 


c. 

+(c,-i-c0 

c._i 


” + >'o 


(AO^X (^-2) 


Contributions from deflection terms 


£j-i £j ^ £j-iC . £jCj_^ 


Cj_,-c^ (c^_i-c,) 


■'j-i 


(B12) 


(OA/tan0o 


Contribution from slope terms 


(i = 1,2,3,...., n) 


Equation (B12) is identical to equation (5c) in the text, and is called the Constant-Domain Nonuniform 
Displacement Transfer Eunction in dual summation form for nonuniform embedded beams based on 
piece wise-linear strain representations. 
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APPENDIX C 


DERIVATIONS OF VARIABLE-DOMAIN IMPROVED DISPLACEMENT 
TRANSFER FUNCTION FOR NONUNIFORM EMBEDDED BEAMS 


Appendix C presents the details of integrations of the beam curvature equation (1) to obtain the 
improved slope and deflection equations for the formulation of variable-domain Improved Displacement 
Transfer Functions for nonuniform embedded beams based on piecewise-nonlinear strain representations. 

Improved Slope Equations 

The slope, tan^(x), of the nonuniform beam in the small domain, < x < x, , between the two 
adjacent strain-sensing stations, {x,_j , x, } , can be obtained by integrating equation (1) as equation (Cl): 


The distributions of the beam depth factor, c(x) [eq. (9)], and bending strain, s{x) [eq. (11)], in the 
small domain, x,_j < x < X; , between the two adjacent strain-sensing stations, {x,_j , x, } , are respectively 
duplicated below as equations (C-2) and (C-3). 


(Based on Piecewise-Nonlinear Strain Representations) 



Slope at X,-. 



(x,._i <X<X,.) 


(C2) 


^(x) = 



(x-x,._i) 


(A/),(A/),,i[(A/),+(A/).,J 



(C3) 


(A/),(A/),,i[(A/),+(A/),,J 


(x,._i <X<X,.) 


Let {A.,6pC, ,^} be defined respectively with equations (C4), (C5), (C6), (C7) as 



(C4) 



(C5) 



(C6) 


^ = (X-X;_i) 


(C7) 
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then the last three terms of equation (Cl) can he written as: 


tan6*(;c) = f 

Jo 




C-i+C,^ 


■d^ + tan6*,_ 


1 ’ 


[0<^<(A/),] 


(C8) 


After carrying out integration of equation (C8), one obtains (ref. 15): 


tan0(x) = 


^log,(c,_, + C,^) + A, - A, ^log,(c,_, + C,^) 




Integration of 1st term in eq. (B-8) Integration of 2nd term in eq. (B-8) 




i(c,_, + C,^)^ -2c,_,(c,_, + C,^) + ctilog,(c,_, + C,^) 


I , C . , C . , 


c, c 


2 ^ 

c;/ 


Integration of 3rd term in eq (B-8) 

A, . 1 fi. 


+ tan0,. 


fi. 


log,(c,_, + C,^) + + -^(c,_, + - 2^c,_,(c,_, + C,a 


c 


-f^-A,%- + B,%llog^c, ,--^cf, + 2^cfi + tan^, , 
I C,. ' Cf ' Cf) ' ' 2Cf ' ' Cf ' ' ' ' 

== -^[cfs._^-A.C.c._,+B,cl^] [log^(c,_| + C,.^)-log^c,_i] + ^^ 


1 B: 


+ 


B 


2 C 


[(C_i + C.^f - cl^ ] - 2 [(c,._i + C,.^) - C;_i ] + tan 0._^ 


" C-i - A<2,C-i + ] [log,(c,._i + C,.^) - log, c,._i ] + 

+ ^^[(C-i + 2c,.-iC,.^ + Cf^)" -cf_i]-2^c,._i[(c,._i + C,.^)-c,._i] + tan6>,_, 

= ^'-1 “ “ log. 

+ (2c, + C,.^^ ) - 2 ^ c,._i^ + tan 


2 Cf 


Cf 


7 ^[cf - A,.C,c,_i + fi,.cf_i ] [log,(c,_i + C,.^) - log, c,._, ] + ^ ^ 

+ ^ (C,-^" - 2c, ._i^) + tan 6>,._, 


(C9) 
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After simplifying, equation (C9) becomes: 


tan 


i + B^cli ]log^ [ 1 + — ^ 


C; 


- 1*1 J 


+ 


1 . 

C, 


A,. + B, 


f n ^ 

i_£Ai 
2 c,. 


(CIO) 


+ tan 


At the strain-sensing station, x^ , we have = x, -x,_j -{Al)- , then equation (CIO) gives the slope 
tan6*,[= tan6*(x,)] at x, as: 


tan0, - A,r,c,^i +B,c,^Jlogjl + -^(A/), 


C. 

(A/), 


C, 


A,+B, 


(A/),- c,,i 

2 C, 


+ tan 0: 


i-l 


(Cll) 


In view of C,- =(c,- — C,_j)/(A/)^ [eq. (C6)], equation (Cll) takes on the form: 


tan 6', = 


(A/)f 


(c,-c,_i)' 


(A/)f 


(c,-c,-i) 


(A/): 

(A/)? 


^ _4 ri thl 

2 


C,-_i 


A, + fi, 


(c-c-.r 


(C-C-i) 

(A/)f 


(A/), 

(A/), c,_j(A/), 


■C-i-A,. 


V 2 c,-c,_iy 
c. - c 


log 
+ tan 0._^ 


+ (AZ),^ 


(A/),- c,_i ; 


A-l 


(AZ), 




(AZ)f 


A, + B,(AZ), 


C - 3c, _i 

2(c,-c,._i)^ 


log,— 

c,._, 


-I- tan 6 : 


(C12) 


After grouping terms, equation (C12) becomes: 

tan e, = 3 [(c,. - c,_i )' f ,_i - A,. ( AZ), c,._i (c,. - c,_i ) + fi, ( AZ)f c," i ] log^ 

(c, 


A 

C._ 


(AZ), 


2(c,-c,_i> 


j[2A,(AZ),(c, - c,_,) + fi,(AZ),^(c, - 3c,_i)] + tan0,_, 


(C13) 


(Z = l,2,3,....,n) 


Equation (C13) is identical to equation (17a) in the text, the variable-domain improved slope equation for 
nonuniform beams in recursive form based on nonlinear strain representations. 
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Improved Deflection Equations 


The deflection, y(x) , can be obtained by integrating the slope equation (CIO) as (ref. 15). 
y(x) = £ tan(x)Jx+ ^ = J^an((^)c?^ + y,_, 


Deflection 

at X: 


= \U 

Jo 


1 ( C 

^[cf - AC,c,_, + ]log, 1 + 

O'; C;_, 


V ^i-1 y 


C; 


A + R 


i_CzL 
2 C; 


+ tan0; 


d^+yi-i 


T[Cf£,_,-A,C,c,_,+B,cL] 


. 0 r, , c, 
c. 


^^l + ^^jlog,^l + ^^j-^ 


1 

H 

c. 


'2 ' I 6 C. 2 j 


+ i^tanR,_j 




/-I 


^ C-I - AC,C-1 + R,cf_i ] 


1 + 


C; 


+il 

A, R, 

ye 

2C, 


1 3 C, j 


c,. >v c,_i 

+ y,.^i + ^tan6»,-_i 


-^^jlog, 1 + ^^j 
/-I J y c-i ^ 


— ^1-^ 


(C14) 


At the strain-sensing station, Xj, we have = x, -x,_j = (A/), , and equation (C14) gives the deflection 
yiV=y{Xi)} as: 


yt = y(^i) = - AC,c,_i + 

^ } 


C, 


-1 


C; 


1 + -^(A/), 


V 


log 


1 + _^(A/), 


V 


-(A/), 


+ 


(A/)? 


2C; 


A,+fi, 


V 3 C, , 


+ y;_i -I- (A/); tan 6>,_i 


(C15) 


In view of C, = (c, -c,_j)/(A/), [eq. (C6)], equation (C15) becomes: 




(A/); 


(C-C-i)' 

m] 


^ -A ^'-' c +Bc^ 

2 ^(-1 / A 7\ ^'-1 ^ 


2(C -C-i) 


(AZ)f 

f 

A +5, 


(A/),- 

(A/), c,_i(A/), ' 


r 1 ^ 


clog. — -(C-C-i) 

C;_, 


3 (C-C-i) 


+ yM+(^)/tan6',.-i 


(C16) 
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After grouping terms, equation (C16) becomes: 


y,- 


(A/)? 




-[(C, - - A,(AZ),(C, - ] 


c,log, — -(C/-c,-i) 


(A/)f 


(C17) 


6(c,-c,_i)- 


{3A(AZ),(c, -c,_,) + B,(AZ)f(c, -4c,_,xJ + y,_, + (A/)Aan0^^^^ 


(Z = l,2,3,....,n) 


Equation (C17) is identical to equation (17b) in tbe text, tbe variable -domain improved deflection 
equation for nonuniform beams in recursive form based on nonlinear strain representations. 

Improved Displacement Transfer Function 

Substitution of slope equation (C13) into deflection equation (C17), and using similar matbematical 
steps of Appendix B, one can write the deflection equation in dual summation form as: 


,4 [iCj - If, - f,-, K-. + B, (A/)’ ] X 

■ ^ 

= Sl ^ c.log,^-(c.-c._i) \-\-yo 

j=i L J 

, (A/), (c , - C ) + fi , (A/)J (c, - 4c , )] 


i-1 

M 


i (AZ), 


k= j+l 


Contributions from deflection terms 
\2 , 


(AZ)^ / (Cj - c,-_i) - AfAl) .c^_fc . - ^ 


;-i 


+ - 


(AZ), 




2 A, ( AZ), (c, - ) + B. (Alf. (c^ - 3c, ) 


Contributions from slope terms 


+ 


I (AO, 


. J=1 


tan6„ 


Contributiohjs from slope terms 


(CIS) 


(Z = l,2,3,....,n) 


Equation (CIS) is identical to equation (17c) in the text, and is called the Variable-Domain Improved 
Displacement Transfer Function in dual summation form for nonuniform embedded beams based on 
piecewise-nonlinear strain representations. 
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Degeneration into Constant-Domain Case 

This section is to prove that when the variahle-domain lengths hecome constant, the variahle-domain 
slope and deflection equations [(C13), (C17)] will degenerate respectively into the constant-domain slope 
and deflection equations [(5a), (5h)]. 


For the constant-domain case, (A/)j - (A /)2 = (A /)3 = = A/, the coefficients, [eq. (C4)] and 6, 

[eq. (C5)] hecome equations (C19) and (C20) respectively: 


^ (A/),,,[2(A/),+(A/)„J.r,_,-[(A/),+(A/),,i]^g,+(A/)f^,,. _ 3g,_. - 4g, + .r,,, 


(A/),(A/),,,[(A/),+(A/).,J 


2AZ 


g - [(A/),. +(A/),.^J^,. +(A/),.^,.^i _ -2s ^ +^,^1 


(A/),(A/),,i[(AZ),+(A/),,J 


2{MY 


(C20) 


Slope Equations 

In view of equations [(C19), (C20)], the slope equation (C13) can he written as: 


tan 9. = ■ 


(A/), 




+ T7-^^^[2A(A/),(c, - c,_,)-t fi,(A/)?(c, - 3c,_,)] + tan0,_, 

^(C; C;_jj 


A. 

C._, 


M 






3£m-4£,+£ 


/-hi 






2£; + 


1+1 1 ^2 
C;_i 


A, (A/) 




log.— 

C,-l 


-H- 


(A/) 


2(c,-c,_J= 
A/ 


f 2)£ — 4£ -I- £ ^ 

JC;_1 ‘tC; -r 


A..(A/) 


(c,-C/-l) + 


^e,_i-2e, + e,,, 


(c,-3c,.-,) 


BiiAlf 


+ tan0; 


2(c,-c,._i) 

M 

"4(c,-c,_,y 

A/ 


2(c, - c,_i )'£,._! + (3 £,._i - 4£,. -t £,.^,)(c, - c,._i)c,._i 
+(£,._! -2£,. + £,.^Jcti 


log.— 

C, 


-1 


-[-2(3£,_i - 4£,. + £„i)(c,. - c,_,) -t (£,_1 - 2e. + £,.,i)(c,. - 3c, ._i)] + tan^,._i 


2(c,-c,_,)^ 
-I- 


(2cf - 4c,.c,_, + 2 c,^i)£,._, + (3£,._i - 4s. + s.^^ )(c,. - c,._i)c,._i 

+(£/-! -2£, + £,-,i)c,^i 

A/ 


log A 

C;-i 


4(c,-c,_,)^ 


[-(6£,_i - 8£, -t 2 £,.^,)(c,. - c,_i) -I- (£,._1 - 2s. + £,+i)(c. - 3c, _i)] -i- tan 6»,_, 
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A/ 




A/ 




A/ 


((2cf - 4c,.c,_i + 2c,^j) + (3c,. - 3c,._| )c._^ + cf_^ ) s._^ 
-{(4c,.-4c;_i)c,._, + 2c,^,)f,. 

(-(6c,.-6c,._,) + (c,.-3c,._|))f,._i 
+ ((8c,. - 8c,._i ) - (2c,. - 6c,._i )) s. + tan 6>,_ 
(-(2c,.-2c,._i) + (c,.-3c,._|))f,.,i 


log,— 


Al 

^2(c,-c,_,)^ 

Al 


(2cf-c,c,_j)c,_i 


(-5 c, -r 3c, j)c, 1 

-(4c,c,_, - 2c,?_i)c, 

1 E AZ 

-t2(3c, -c,_i)c. 

_ + C,C,_,C,^j 

c,-_i 4(c,-c,_i) 

_(-c,--c,_i)c,,i _ 


+ tan 


[(2c, - c,._i)c,f,_i - 2(2c,. - c,_i)c,._i^, + c,c,._i^?,^i]log, ^ 

^i-l 

-[-(5c, - 3c,._i)f,_i + 2(3c, - c,_i)^?, - (c, + c,_i)f,^i] -f tan^,.. 


After grouping terms, equation (C21) becomes: 

Al c 

= t; ^[(2c, -c,_,)(c,f,_i -2c,_i^,) + c,c,_i^r,^i]log,^ 

2(c,-c,_i) c,_i 

^^-^[(5c, - 3c,_i)f,_, -2(3c, -c,_i)^r, +(c, +c,_i)iT,^i] + tan^,_i 
4(c,-c,_i) 


Equation (C22) is identical to the improved slope equation (7a) for the constant-domain case. 
Deflection Equations 

In view of equations [(C19), (C20)], the deflection equation (C17) can he written as: 


y, = 


+ - 


(AZ)f 

(Al)^ 


- [(c,- - c,_i f c,_i - A, (Al). (c, - c,_, )c,_i -f B, (A/)," c,?_i ] 


c,log, — -(c,-c,_i) 


6(c,-c,_i) 


y [3A( AZ), (c, - c,_i ) + B, (Al)^ (c, - 4c, ), ] + + (Al), tan B,_, 


(C21) 


(C22) 
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X 


{Mf 






3e,_i-4£, + £,, 




£,_j - 2£,. + £,.+1 I 2 

^i-1 


A,- (A/) 


5,(A/r 


■ 

( C ^ 


C;l0g, 

S 

-(C;-C,_i) 


V ‘'i-l 7 



+ ■ 


(A/)^ 


2(c,-c,_i) 


(A/)^ 


3£--i 4£. + ^ f £,-_i 2£. + A c. 4c-_^ 

2 j L 2 )Xc.-c,_,) 


A, (AO 
\2 


a (AO" 


+ + (A/),. tane,._j 


2(c, - c,_i)' + (3f;_i - 4f . + £,^i)c,-_i(c, - c,_i ) 


2(c, C;_j) - 2s- + £,-+i)cf_j 






+ - 


{Mf 


I2{c,-c.,j 


■-3(3f,_i-4f, + f,-^i)(c,-c,_i)' 
_+(f,_i-2f, + f,^i)(c,-4c,_i) ^ 


+ 3',-i + (A/),tan6',_i 


{Mf 


2{c-c,J 


+ 


(-4(c,-c,_i)c,_i-2c,^i)£, 

c,log, — -(c,-c,_i) 
*^(-1 

((C,-C,_i)c,_i+C,"_l)£,,l 



{Mf 


{Mf 


2(c,-c,_.)^ 

{Mf 


(-9(c, c, + (c,.- 4c, ._! ))£,_! 

+ (1 2(c,. - c,._, ) - 2(c,. - 4c,._, )) £,. 
(-3(c,. - c,. _i ) + (c,.- 4c,. _i ))£,.,! 

[{2cf - c,c,_i )c,_i - 2(2c,c,._i - cf_i )c, + c,-c,_iC,^i ] 


12(c,-c,_j' 


+ j,._i+(A/),.tan6>,._i 


c,log,^-(c,-c,_i) 
c,- , 


12(c,-c,_jf 
After grouping terms, equation (C23) becomes: 
{Mf 


[-(8c, - 5c, _i )c,_i + 2(5c, - 2c,_i )c, - (2c, + )c,^i ] + y,., + (A/), tan 


yi = 


2(c,-c,_,y 

{Mf 


-[(2c, - c,_i)(c,c,_i - 2c,_iC,) -t c,c,_,f,^j] 


c,log, — -(c,-c,_i) 
c,--i 


— — ^[(8c, - 5c,_i)c,_, - 2(5c, - 2c,_i)c, + (2c, + c,_, )c,^, ] + y,_, + (AZ), tan 

iztc, c,_,j 

Equation (C24) is identical to the improved deflection equation (7b) for the constant-domain case. 


(C23) 


(C24) 
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Displacement Transfer Function 

Combining equations (C22), (C24), and after grouping terms, one obtains the deflection equation in 
dual summation form as: 


J=l 

2{C.-C.,f ^ ^ ] X 

Cj 

X c.log,- — (c-c i) 

L ‘ \ 

12(c.-c. 2c._i)£.+(2c.+c._i)e.,i] 

>+yo 

7=1 

Contribution from deflection terms 

3 [( 2 c/ 2c._,£.) + c.c._ie.,i]log^ ^ 

Cj-i) 

4(c.-c. 2(3 c c,_i)e^.+(c, + c,_Je^,i] 

(C25) 


+ (i)Al tan0Q 

Contributions from slope terms 

Equation (C25) is identical to equation (7c), the Constant-Domain Improved Displacement Transfer 
Function in dual summation form for nonuniform beams based on piecewise-nonlinear strain 
representations. 
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APPENDIX D 

DERIVATIONS OF THE VARIABLE-DOMAIN LOG-EXPANDED 
DISPLACEMENT TRANSFER FUNCTION FOR SLIGHTLY 
NONUNIFORM EMBEDDED BEAMS 

(Based on Piece wise-Nonlinear Strain Representations) 

Appendix D presents the mathematical derivations of the slope equation (18a), deflection equation 
(18h), and Displacement Transfer Function (18c) for slightly nonuniform beams through expansions of 
logarithmic functions. 


Log-Expanded Slope Equations 

The slope equation (17a) is duplicated helow as equation (Dl). 


tan 9: = 


(A/), 


■ [(c, - c,._i f - A. (A/), c._, (c. - c,._i ) + B, (Al)f cl, ] log^ 


(c,-c,-i)- 

+ T7-^^^[2A(A/),(c, - c,_,) + fi,(A/)f(c, - 3c,_,)] + tan0,_. 


(Dl) 


For slowly changing c(x) [that is, (C(/C(_i) 1], the logarithm term, log^(c, /c,_j) , in the slope equation 

(Dl) can he expanded in series form up to the third order term as (ref. 15): 

rc„ A irc,. Uc,. 


1 I J ^ I J ^ I 1 I 


C, - c,_ 


C:_ 


I C/-g,-l I (g, -c,-i) 


6cr 


2c, -I 3cf_i 

-\_6cl, - 3c,_i(c,. - c,._i) + 2(q. - c,_,f - ....] 


6c, -I 


i[6cA -(c,. -c,._i)(3c,_i -2(c,. -c,_i)) -....] 




6c 


[6c,^i + (C; -c,._i)(2c,. - 5c, ._! )-....] 


(D2) 
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In view of equation (D2), equation (Dl) can be written as: 


tan 9 : 


(A/), 


- [(c,- - c,_i f - A, ( AZ), c,._i (c, - c,_i ) + B, ( AZ)f c]_, ] x 


X % + (c,- - c,-i)(2c,- - 5c,_i)- ....] 

6c, -1 


+ 9 / ' x 2 [2A,(AZ),(c, c,_,) + B,(AZ)f(c, 3c,_i)] + tan6',_i 

^vC, c,_,; 


(c, - c,_,)"[6c,"_i + (c, - c,_i)(2c, - 5c,_,)- ....]f,_, 

(AZ), 

-A,(AZ),c,_i(c, - c,_i)[6c,"_i -t (c, - c,_i)(2c, - Sc,.,) - ....] > 

6c,^,(c,-c,_,)^ 

+B,(AZ)f c,^i [6c,^i -t (c, - c,_i)(2c, - Sc,.,) - ....] 

, (A/), 

^[6A,(AZ),c,^,(c, - c,_,) + 3B,.(AZ),^c,^.(c, - 3c,_,)]+ tanB,_ 


6c,._,(c,-c,_i) 


(AZ), 

6c,^i(c,-c,_i)" 


(AZ), 


6c,^,(c,-c,_,)^ 


(AZ), 


6c,_,(c,-c,_,) 


(c, - c,_,)" [6c,"_i + (c, - c,_i)(2c, - 5c,_i) - 
-A,(AZ),c,_j(c, - c,_i)[6c,^i + (c, - c,_,)(2c, - 5c,_,) - ....] 
+B,(AZ)fcL[6ct, + (c, -c,_,)(2c, -5c, 

[6A,(AZ),c,^,(c, - c,_,) + 3B,{Mf,cUc, - 3c,_,)] 

(c, - c,_,)" [6c, ", + (c, - c,_,)(2c, - 5c,_,)- 
-6A,(AZ),c,^,(c,-c,_,)-A,(AZ),c,_,(c,-c,_,)^(2c,-5c,_,)-..^ 
+6B,(AZ)f c,^ + - c,_,)(2c, - 2c,_,) - 3c,_,) - .... 

6A,(AZ),c,^,(c, - c,_,) + 3B,{Mtcl ((c, - c,_,) - 2c,_,) 

(c, - c,_i)" [6c,"_i + (c, - c,_i)(2c, - 5c,_, )....]f,_i 
-A, ( AZ ), c,_i (c, - C;_, )^ (2c, - 5 c,_, ). . . . [ + tan B,,, 

+2B,(AZ)^c,^,(c,-c,_,)l... 


+ tan B,_, 


+ tan 9 : 


(D3) 


After cancelling (c, terms, equation (D3) becomes: 


tan 0, = 


(AZ), 


(6c,"_i + (c, - c,_,)(2c, - 5c,_,))f,_i 
-A,(AZ),c,_,(2c, -5c,_,) + 2B,(AZ)fc,^, 


+ tan 0,_j 


Equation (D4) is tbe Log-expanded slope equation (18a) in the text. 


(D4) 
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For uniform beam = C, = C), equation (D4) becomes: 


tan 6*. = 


(A/), 


A,. 


B, 


s,_,+^(AlX+^(Al)t 


+ tan 6',_i 


which is identical to equation (E8) in Appendix E. 


(D5) 


Log-Expanded Deflection Equations 


The deflection equation (17b) is duplicated below as equation (D6). 




(AZ) 




- [(c, - c,_i f - A. ( AZ), (c, - c,_i )c,_i + B. (Al)f cf_i ] 


c,log, — -(c,-c,_i) 


(AZ) 


6(c,-c,._,)- 


- [ 3A( AZ), (c, - c,_, ) + B, {Alt ic. ~ 4c, ), ] + + (AZ), tan d,_, 


(D6) 


Eor slowly changing c{x) [that is, (c,/C(_i) ^ 1], the logarithm term, log„(c, Zc,_,) , in the deflection 
equation (18b) can be expanded in series form up to the fourth order term as equation (D7) (ref. 15): 


, c, _rc, if c, y if c, ,v ifc, Y 

°^'c,_i U,_, J 2U,_1 J ^3U,_1 J dU,., ) 


I2ct, 


-[l2c,^i - 6c,"_i(c, - c,_i) + 4c,_,(c, - c,_i)" - 3((c, - c,.,)' + ....] 


Then the term, [c,- log,(c,- / C,_j) — (c,- — C,_j)] , in equation (D6) can then be written 


as: 


(D7) 


c, 


Cilog,^ — (c,-c,-i) 


= 4^''~‘^ [12c,^_i - 6c,"_i(c, - c,_i) + 4c,_i(c, - c,_i)" - 3((c, - c,_i)^ + (c, - c,.,) 

12c,_, 


c,(c,-c,_i) 


l2cU 

C,(c, - c,_i) 


l2cU 


1 2r 

12c, - 6c,^i(c, - c,_i) + 4c,_i(c, - c,_,)" - 3((c, - c,_,)' ^ 

. 

1 2c^ 

— (c, - c,-i ) - 6c,", (c, - c,_i ) + 4c, (c, - c,_i )" - 3(c, - c,_, )" 






/-I 


Y— ^^[6c,^i +(4c,c,_i - 6c,^i)(c, -c,_,)- 3c,(c, -c,_,)"] 


(D8) 
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In view of equation (D8), equation (D6) can be rewritten as: 




(A/)f 






I2cl 


[6c,^i + c,._i(4c,. - 6c,._i)(c,. - c,_i) - 3c,. (c,. - c._^f ] 


(A/)f 


_ ^[3A(A/),(c, - + - 4c,_,)] + + (A/), tan^,_, 

Otc, C,_jl 

(c,. - c,._i f [6c,^i + c,._i(4c,. - 6c,._, )(c,. - c,._i) - 3c,. (c,. - 
-A.(A/),.c,._i(c,. - c,._,)[6c,:Li + c,._i(4c,. - 6c,._,)(c,. - C;_,) - 3c. (c, - c,._ 
+B,(A/)f c,"_i \_6cli + c,._i(4c,. - 6c,._i)(c,. - c,._i) - 3c,. (c,. - c,._i)"] 

r[6A,.(A/),.c,:‘_i + 2B,.(A/)f c,^‘_i(C; - 3c,._i )] + + (A/),, tan 


(A/)f 

12c,^,(c,-c,,)^ 


(A/)f 


(A/); 


12 c,^,(c,-c,._J^ 

(c,. - C ;_1 f [ 6 c,^i + c,._i (4c,. - 6 C ;_1 )(c,. - C,._, ) - 3c,. (c, - c,._i f ] 
-A.(MX(c, - c,._i)c,._i [ 6 c,i^i + c,._i(4c,. - 6 c,._i)(C; - c,._i)- 3c,.(c,. - c._, 
+6A,(MXcUc,-c,J 

+B,(M). cf_^ [6cl^ + c,._;(4c,. - 6c,._i)(c,. - c,_,)~ 3c,. (c,. - 
+2B,(A/),^c,^,[(c,.-c,_,)-3c,_J 

+ y/-i + (A/),tan^;_i 


12c,^,(c,-c„) 


(A/)^ 


12c,:L:(c,-c,_,)^ 


(c,. - C;_i f [6c, + c,._i (4c,. - 6c;_i )(c,. - c,._, ) - 3c,. (c,. - c,._i f ] 

6c,^!., + c,._.(4c,. - 6c,._.)(C; - c,._. )- 3c,.(c,. - c,._. 

-A,(A/),(c,-c,_Jc,_, 


-6c: 


6cl^ + c,._| (4c,. - 4c,._i - 2c,._i )(c, - c,._i ) - 3c,. (c,. - C;_i ) 
_+2c,"_,(c,.-c,._i)-6c,^, 

+ y,._, + (A/),.tan6';_i 

(C; - c,_i)" [6c,?_i + c,.^i(4c,. - 6c,.^i)(c,. - c,._i) - 3c,. (c,. - c,._i)" 

4 , ;ri-A,(A/),.(c,.-c,._i)c,..i[c,^i(4c,.-6c,^i)(c,-c,_i)-3c,.(c,.-c,._if] 

12c,._i(c,.-c,..i) ^ 

+fi,(A/)fc,^,[4c,„.(c, -c,.,f -3c,(c, -c,_,)^] 


(A/)? 


_ (A/),^ 
12c,^ 


+ y,.^i + (A/),.tan0,._; 

[6c,^L, + c,._i (4c,. - 6c;_, )(c,. - c,_i ) - 3c; (c,. - C;_, f ] f,._, 
-A,.(A/),.C;_i [c,._i(4c,. - 6c,._i) - 3c,. (c,. - c,_i)] 
+fi,(A/)?cL(4c,_,-3c,) 


^ + }',_i + (A/),.tan6',._i 
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(A/)f 

12cf, 


(A/)f 

12c^, 


(A/)f 

12cf, 


[6cti + c,._i(4c,. - 6c,_i)(c,. - c,_i) - 3c,. (c, - 
-A;(A/),.c,_i[c,._i(4c; -4c,_i -2c,._i)-3c,.(c; -C;_i)] 
+B,(A/)^c,:i,(4c,_j-3c,) 

[6c, + c,._i(4c,. - 6c,_i)(c,. - c,_i) - 3c, (c, - c,_i)"]c,_, 
-A,(A/),c,_, [4c,_,(c, - c,_,) - 2c,:ii - 3c, (c, - c,_,)] 
+fi,(A/)fc,^,(4c,.,-3c,) 

L 

[6c,li + c,_i(4c, - 6c,_i)(c, - c,_i) - 3c, (c, - c,_,)" ]e,_, 
+A,(A/),c,_,[2c,^, + (3c, -4c,_,)(c, -c,_,)] 
+B,(A/)fc,^,(4c,_,-3c,) 


^ + y,_i+(A/),tan6',_i 


■ + yM + (A/),tan6>,_i 


^ + y,-i + (AO,tan0,_i 


(D9) 


Term L in equation (D9) can be rewritten as: 

L = [6c,^i +c,_i(4c, - 6c,_i)(c, -c,_,)- 3c,(c, -c,_if] 

= [6c,^i + (4c,c,_, - 6c,^i - 3c, (c, - c,_,)}(c, - c,_, )] 

= [6c,^i + (4c,c,_i - 6c, - 3c," + 3c,c,_i )) (c, - c,., )] 

= [6c,^j + (7c,c,_, - 6c, - 3c,")(c, - c,_,)] 

= [6c,^i + (7c,"c,_i - 6c,c,"_i - 3c," - 7c, c,^, + 6c, 1, + 3c, "c,., )] 

= [6ct, + (I0cfc,_, - 13c, ct, - 3cf + 3ct, + 3ct,}] 

= [6c,^i + (lOc,"c,_i - 13c,c,"_i + 3c,"_i - 3(c," - c,".,)}] 

= [6c,^i + (lOc,"c,_i - 10c,c,"_j - 3c,c,^i + 3c, - 3(c," - c,"_,))] 

=[6c,"i +{lOc,c,.,(c, -c,.i)- 3c,"_i(c, -c,.,)- 3(c," -c,"i))] 

= [6c,?, +c,,(10c, - 3c,,)(c, -c,,)- 3(c," -c,^,)] (DIO) 
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In light of equation (DIO), equation (D9) takes on the final form as: 


Ji = 


(A/)f 

12cf, 


[6cti + c,._i(10c,. - 3c,._i)(c,. - - 3(cf - 
+A,.(A/),.c,.^i [2c, + (3c,. - 4c,._i)(c,. - c,.^i )] 
+5,(A/)?c,^,[c,_,-3(c,-c,^,)] 


> + y,_i + (A/),.tan6>,._i 


(Dll) 


Equation (D1 1) is the Log-expanded deflection equation (18b) in the text. 
For uniform beam (C,_j = C, = C), equation (Dll) becomes equation (D12): 


(A/)^r A. B , 

^,_, + ^(A/),+^(A/),^ 
2c 3 6 


+ y,._i + (A/),.tan6',._ 


which is identical to equation (El 1) in Appendix E. 


(D12) 


Log-Expanded Displacement Transfer Function 

Combining equations (D4) and (Dll) and using similar mathematical steps to Appendix B, one can write: 


1 ^ (A/)^ 
= 12^ c^ ^ 

,;=I ^i-\ 

+A,(AZ),[2( 

10c . - 3c,._, )(c . - c ) - 3{c] - 
’■_i+c,._i(c.-c._i)(3c.-4c._i) 
[c._,-3(c.-c,_,)] 

A. 

^ + >'0 

c 

;=1 [_A:=J+1 J 

!!ontributions from deflection terms 

6c^i + (c^.-c._i)(2c^.-5c^._j))e^._i 

-A^. ( AZ) . c,_, (2c, - 5 c^._, ) + 2B. {M)] c]_, _ 

r n 1 (D13) 

+ X(AZ)j tan 00 i 


Contributions from slope terms 


Equation (D13) is equation (18c) in the text, the Variable-Domain Log-expanded Displacement Transfer 
Functions in dual summation form for slightly nonuniform (including uniform) embedded beams based on 
piecewise-nonlinear strain representations. 

Degeneration into Constant Domain Cases 

For constant-domain length, (A/)j = (A /)2 = (A /)3 == == A/ , coefficients A, [eq. (C4)] and B- 

[(eq. (C5)] become equations (D14) and (D15): 


^ ^ (A/)„,[2(A /), + (A/),,J^,_, - [(AZ), + {M),Js, + _ 3s, - As, + 


(AZ),(AZ),,,[(A/),+(AZ)„J 


2AZ 


(D14) 


S = 


- [(A/),. + (AZ).^Jf,. + (A/),, S;_, - Is. + 


(AZ),(AZ),,,[(AZ), + (AZ), 


/+! J 


2(A/)' 


(D15) 
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Slope Equation [(A/)^ A/] 


In view of equations (D14) and (D15), slope equation (D4) can be written as: 


tan 9. 


(A/), 


12c^ 


(A/), 


+ 2(c,. - c,._i)(2c,. - 


+6A,<,(A/),-4A,(A/),c,_,(c,-c,_i) + 45,4,(A/)f 


+ tan0. 


/-I 




6 cL 


(2c,.-5c,._,)^r,._i 


+ ^(A/), - - c,_,) + ^(A/)f 


+ tan 9- 


Ci-x 




6c?, 


(2c,.-2c,._i-3c,._i)^r,._i 


3^r,._i - 4^r,. + ^ 3s. - As, + , s,_, - 2s, + s,^, 


6c, 


-(c,-C,_,) + - 


-1 


+ tan 9, 


AL 

c,--i 


s,_i + 


2(c,-c,._,) 3(c,-c,._i) 


6cf_i 


6c, .X 


3s,_, - As, + s,,, ^ 3s,_,-As,+s„, ^^ _ 3s,_,-As, + s„. 


6c, -1 


+ 


s,_,-2s, + s,^, 


+ tan 9,_, 


AL 

C /-1 


A/ 


^1-1 


M 


C:_ 


^i-X + 


(.c,-c,_,y c, lie, . . 

— ' — s, I —S, , +-S, I S, I + —s, I S, I + — I 

3cf_, 2c, 2 '■* 4 '■* 2c, _i 2 6 


1 


1 


-As, -As, -As, -2 s, 

L + L(. L + L 

4 6c, ' 6 6 


£i+x_ , £j+x _ £j+x , £j+x 


4 6c, _i ' 6 

r, (c,-c, ,)' 3 1^ f, 2 I 2c, ^ ( I c, ^ 

J I 4 6c,_J 


+ tan 9,_, 




3cti 4 6; '"^ 3 3 3c,_,; ' I 4 6c, _ 


I 

7 5 4(c,-c,_if^ T4 2c, ^ f 3 2c,) 


+ tan 9,_, 


12 ^ 12c, 


A J^‘'‘^l3 3c,_J^'^l 12 ^ 12c,_J 


+ tan 9, 


1-1 


(D16) 


After grouping terms, equation (D16) becomes: 


tan 9, = -^[[5cA + 4(c, - c,_i f]s,_, - 8c,_, (c, - 2c,_, ) + c,_, (2c, - 3c, )s,^, } + tan 9,_, 
lzc,_[ 


(D17) 


Equation (D17) is identical to equation (8a) in tbe text, and thus confirming the mathematical accuracy of 
the Log-expanded slope equation (D4). 
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Deflection Equation [(A/); A/] 


In view of equations (D14) and (D15), deflection equation (Dll) can be written as: 


y,' = 


(A/)f 

\2cU 


12ct: 


'[6ct, +(10c,c,^, - 3c,^,)(c, -c,,)- 3(cf 
+A,.(A/),.c,_i [2c, + (3c,. - 4c,.^i)(c,. - c,..;)] 
+fi,(A/)fc,^,[c,.,-3(c,-c,.,)] 

+ y,._, +(A/),.tan6>,._i 

[6c,i^i + (10c,.c,._i - 3c," 1 )(c,. - c,._i ) - 3(c," - cl ^ )] 

c._i [2c,"_ 1 + (3c,. - 4c,._i)(c,. - c,._i)] 


+ 


€._^- 2 s ,+ s ,_ 


-]c,^i[c,._i-3(c,.-c,._i)] 


> + j,._i+(A/)tan6';_i 


(A/)" 

12c, ^ 


> + y,-i+(A/)tan6»;_i 


(A/)" 
' 24cf 


(A/)" 

24c,^ 


[6c,"_i + (10c,c,._i - 3c,"_ i)(c, - c,_i) - 3(cf - c,"_i)]^,_i 
c,._i [2c,"_i + (3c,. - 4c,_i)(c, - c,_i)] 

c,._i [2c,"_i + (3c, - 4c,_i)(c, - c,_i)] 

c,-i [2c,"_ 1 + (3c, - 4c,_i)(c, - c,_i)] 

c,"_i [c,_i - 3(c, - c,_i)] 

c,"_i [c,_i - 3(c, - c,_i)] 

+ [^]c,"_i[c,_i-3(c,-c,_i)] 

[l2c,"_i +(20c,c,_i -6c,"_i)(c, -c,_i)-6(c," -c,^i)]^,._i 
-3^,_iC,_i [2c,^i + (3c, - 4c, _i )(c, - c,_i )] + [c,_, - 3(c, - c,_i )] 

+4^,c,_i [2c,"_i + (3c, - 4c,._i)(c, - c,_i)] - 2^,c,"_, [c,_i - 3(c, - c,_,)] 
[+2c,^i + (3c, - 4c,_i)(c, - c,_i)] + [c,_i - 3(c, - c,_,)] 

[l2c,"_i +c,_i(20c, -6c,._i)(c, -c,._i)-6(c," -c,"_i) 

-6c,"_i - c,_i(9c; - 12c,._i)(c, - c,_i) + c,"_i - 3c,"_i(c, - c,_i)]c,_i 
+c,_i [ 8c,"_, + (12c, - 1 6c,_j )(c, - c,_i ) - 2c,"_i + 6c,_, (c, - c,_, )] c, 

-c,_i [2c,"_i + (3c, - 4c, )(c, - c,_i) - c,"_i + 3c,_i(c, - c,_, )]c,^i 


+ +(AZ)tan6',_, 


+ y,-i + (A/)tan6’,_i 
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{Mf 

2AcU 


\jcl, + C._, (1 Ic. + 3c,._i )(c,. - c,._i ) - 6 (cf - cf _, )] £;_1 

Fi 

+ C;_i[ 6 c,^i +( 12 c, - 10 c,_j)(c, -c, 

^2 

[cf_i + (3c,. - C ;_1 )(c,. - c,._i )] £,.^1 


+ j,._i+(A/)tan0,._j 


(D18) 


Ther terms {Fi, F 2 , F 3 } in equation (D18) can be rewritten as: 


F, - [7ct, + c,_,(l Ic, + 3c,_J(c, - c,_.) - 6(cf - 

= [7c;^. + 1 lcfc,_, - 1 lc,c;i, + 3c, cl - 3cl - 3cf + 3cl)~ 3(cf - c^j] 

- [icl + 8cfc,_, + 3cfc,_, - 8c, c^ - 3cf - 3(cf - c;^,)] 

- [7c,:^i + 8c,c,_,(c, - c,_,)- 3c,"(c, - c,_,)- 3(c,^ - c,^,)] 

= [7c, + (c, - c,_, )(8c,c,_. - 3c,^ ) - 3(c,^ - c,^, )] 

F 2 = c,_i [6c,^i + (12c, - 10c,_,)(c, - c,_i)] 

- 2c,_i [3 c,!1i + 6cf - 5c,c,_| - 6c,c,_, + 5c,"_, ] 

= 2c,_i[8c,!i, + 6c,"-11c,c,_i] 

= 2c,_i [8c,"_ , + 3c," + 3c," - 8c,c,_, - 3c,c,_, ] 

= 2c,_i [3c," + (3c, - 8c, )(c, - c,_,)] 


^3 = 


-c,_i [c,"_i + (3c, - c,_i)(c, - c,_i)] 

-c,_i [c,^i + (3c," - c,c,_, - 3c,c,_i + cl)~\ 
-c,_i [c," + (2c," - 4c,c,_, + 2c,"_,)] 


= -c. 


i-l 


c-+2{c-c,_,f] 


(D19) 


(D20) 


(D21) 


Substitutions of equations (D19), (D20) (D21) into equation (D18) yields: 

_ (A/)" |[7c,"_i + (c, -c,_,X8c,c,_,-3c,")-3(c,"-c,"_,)]c,_i 
24 c,"i |+2c,_,[(3c," +(c,-c,_i)(3c,-8c,_,)]c,-c,_i[c," + 2(c,-c,_i)"] 


+ y-_^ + A/ tan^,_. 


(D22) 
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Equation (D22) is identical to equation (8b) in the text, and thus confirming the mathematical accuracy of 
the Log expanded deflection equation (Dll). 


Displacement Transfer Function [(A/); A/] 

Combining equations (D17) and (D22), and after grouping terms, one obtains the displacement transfer 
function in dual summation form as: 




(Alf 


24 -c._. 


[lc]_, + (c, - c^_A(8cc_, - 3cp- 3(c] - c]_,)\e^_, 
+2c [(3c^^ + (c . - c._^ )(3c^ - 8c )] £, 




f+yo 


Contributions from deflection terms 


(D23) 


12 i Uc,_,(2c,-2c,.,)e,„ < 


Contributions from slope terms 


Equation (D23) is identical to equation (8c) in the text, the Constant-Domain Log-expanded Displacement 
Transfer Eunctions in dual summation form for slightly nonuniform (including uniform) embedded beams 
based on piecewise-nonlinear strain representations. 
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APPENDIX E 

DERIVATIONS OE VARIABLE-DOMAIN IMPROVED DISPLACE 
TRANSEER EUNCTIONS EOR UNIEORM BEAMS 


(Based on Piecewise-Nonlinear Strain Representations) 


Appendix E presents the derivations of the slope and deflection equations for the special case of the 
uniform beam based on piecewise-nonlinear strain representations. 

Slope Equations [c(a:) = c] 

The slope, tan^.r) , of the uniform beam [c(.r) = c] in the small domain, x,_i < x < x, , between the 
two adjacent strain-sensing stations, {x,_j , x, } , can be obtained by integrating both sides of equation (El): 


The piecewise-nonlinear equation (16) representing the bending strain, £’(x) , within the small domain, 
X;_i < X < X, , between the two adjacent strain-sensing stations, {x,_i , x, } , is duplicated below as equations 
(E2) -(E5): 



Integration of right hand 
side of equation (1) 



Slope at Xj_^ 


Intergation of left hand side of equation(1 ) 


Si.^) = Si-i 



(A/),(A/),,i[(A/),+(A/),,J 



(E2) 


(A/),(A/),,i[(A/),+(A/),,J 


(x,._i <X<X,.) 


Eet, 



(E3) 



(E4) 


<^ = (x-x,._i) 


(E5) 
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In view of equation (E2), (E3), (E4), and (E5), the last three terms of equation (El) can he written as: 


tan0(^) = -f(f,^i+A,^ + B,^')fii^ + tan0,^i ; [0<^<(A/),] 

c •'0 


(E6) 


After carrying out integration of equation (E6), one obtains (ref. 15): 


tan^(x) = — 
c 




+ tan^,_i ; [0<^<(A/),] 


(E7) 


At the strain-sensing station, x, , we have i^==X;-x,_j then equation (E7) gives the slope, 

tan6*,[= tan^(X;)] at x, as: 


tan^, = 


— +^(V>n + tan(^._, 
c 2 3 


(/ = l,2,3,....,n) (E8) 


Equation (E8) is the slope equation in recursion form. 

Deflection Equation [c(x) = c] 

Deflection can he obtained by integrating the slope equation (E7) as: 


y(x) = r tan6>(^)<i^ -t y,_i = -f< 

JO (f Jo 


e t: 2 e 3 


tan^,_i|<i^-ty,_i (E9) 

[0<^<(A/),1 


Carrying out integration of equation (E9), one obtains: 


2c 


A- B- 2 

3 D 


y,._i+^tan^,_i ; [0<^<(A/),] (ElO) 


At the strain-sensing station, x, , we have = x, -x,_j = (Al)- , and the deflection, y, [- ^ is 

given by equation (ElO) in the form: 


(A/)^ 

2c 




-1 


-(A/), 


A(A0? 

D 


-t(A/), tan6»,_i 


(t=l,2,3,....,n) (Ell) 


Equation (Ell) is the deflection equation in recursion form. 

Displacement Transfer Function [c(x) = c] 

Substitution of slope equation (E8) into deflection equation (Ell), and following the mathematical 
steps presented in Appendix B, one obtains the deflection equation written in dual summation form as: 
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(E12) 






Contributbn from deflectionterms 


;-l 

:Z- 


Z<A'). 

k=j+l 


{Aiy 


A/ 5, 2 

^2_,+^(A/)2+^(A/)5 


Jo 


Z<")j 

j=l 


i&nOr, 


Contributbns from slope terms 


=0 for cantilever beams 


Equation (El 2) is called the variable-domain Improved Displacement Transfer Eunction in dual 
summation form for uniform beams based on piecewise-nonlinear strain representations. 

Degeneration into Constant-Domain Cases 

Eor constant-domain length, (A/)j = (A /)2 = (A /)3 = = A/ , coefficients [eq. (E3)] and B- 

[(eq. (E4)] become equations (E13) and (E14): 


A, = 


(A/),.,i[2(A/); +(A/),,J^r,_i - [(A/),. +(A/),,j2.t,. +(A/)f^r,,i 3^r,_, -4^r, +^r. 


(A/),.(AZ),,,[(AZ),+(AZ),,,] 


2AZ 


(E13) 


g - [(A/),. +(A/),.^J^,. +(AZ ),.^,.^1 _ -2s, 


(A/),(A/),,i[(AZ),+(AZ),,J 
Slope equation [c(x) = c] : 


2{MY 


In view of equations (E13) and (E14), slope equation (E8) can be written as: 

(AZ), 


tan 9: = ■ 


^/-i+y(AZ),+^(AZ)f 


-I- tan 9,_, 


c 


+ 


( ^s,^^-4s,+s,^l^ 

AZ 

s,^, -2s, +s,^, 

[ 2AZ ; 

2 

[ 2{Alf ) 


(Aiy 


-f tan 9,_, 


= ^ - 3(3f - 4s, + s,^i ) + 2(e,^, - 2s, + ^ )] + tan 9,_, 
12c 

= -9^,-1 +12^, -3s,^,+2s,_, -4s, +2s,^,) + tan9,^, 

12c 


(E14) 


(E15) 


After grouping terms, equation (E15) becomes: 


tan 9, = — (5f ,._i -h 8^ , -•&,+!) + tan 9,^, ; (Z = 1,2,3, 

12c 


(E16) 


Equation (El 6) is the degenerated form of equation (8a) in the text for C,_, = C, = C. 
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Deflection Equation [c(x) = c] 

In view of equations (E13) and (E14), deflection equation (Ell) becomes: 


yi 


mf 


2c 


(Aiy 


3 6 


+ +(A/),. tan6’,._i 


2c 


(Air 

24c 
_ (Alf 

24c 


y-i + 


f 3c,._i -4c,. 

Al 

C,.^1 -2c,. +c,.^i 

1 2A/ ; 

3 

1 2{Alf ) 


(Air 


+ y + A/ tan 0. 




[l2f - 2(3f - 4e. + ) + (£.^, - 2£. + )] + y.^, + Al tan 0._^ 

[l2f,_i - + 8f,. - 2f,.^i + £,_i - 2s, + €,^,]+ y,_, + Al tan (9,_, 


After grouping terms, equation (E17) becomes: 


(E17) 


24c 


(El 8) 


Equation (El 8) is the degenerated form of equation (8b) in the text for = C, = C. 

Displacement Transfer Function [c(x) = c] 

Combining equations (E16) and (E18) into one equation, and after grouping terms, one obtains: 

(Air ' (Air ‘ ' 

X- = + “^;+i) + )'o + 8e,- -£;+i)] + (OA/tan0o 

Contributions from deflection terms Contributions iiom slope terms 

(/ = l,2,3,....,n) 


Equation (E19) is the degenerated form of equation (8c) in the text for 0,^^ = C, = C. 

Special Case 

Eor a special case when the strain distribution along the uniform beam is linear, one can write equation 
(E20), 

(E20) 

or equation (E21), 

^,+1 -c,_i (E21) 
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In light of equation (E21), the slope equation (El 6) becomes: 


tan 6: = -^(5s , , + Sf,. - ) + tan 0. , 

l 1 ^~1 I H"l ^~l 

I2c 


A/ 

I2c 


[5 + 8f,. - (2f; - )] + tan 0._^ 


After grouping terms, equation (E22) becomes: 


A/ 


tan = — (f +£i) + tan ; (t = 1,2,3, ... . ,n) 

2c 


Equation (E23) is tbe degenerated form of equation (6a) for C■_^ = C- = C. 
Also, in view of equation (E21), tbe deflection equation (E18) becomes: 


24c 


\lSi_i + 6s I - (2s^ - Si_i )] + + AZ tan d-_^ 


After grouping terms, equation (E24) becomes: 


(Mf 


Ji = —— (2^,_i +£i) + y,._i + A/ tan6Z,._i ; (i = 1,2,3, . . . . ,n) 

DC 


Equation (E25) is tbe degenerated form of equation (6b) for = Cj = C. 
Einahy, in view of equation (E21), tbe deflection equation (El 9) becomes: 
(Alf ^ ^ ^ (Alf 


yi - 24 ^ + + “^;+i)] + (OA/tan0o 


(AZ) 


24c 
(Alf 


Contributions from deflection terms 
2 i 


Contributions from slope terms 


X [7e,._, + 6c . - (2e . - e )] + 


+ ^ [(Z - j3[5e ._, + 8e . - (2e . - e ._,)]] + (Z)AZ tan ©0 

12c 


(E22) 


(E23) 


(E24) 


(E25) 


(E26) 
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After grouping terms, equation (E26) becomes: 


y. - ^ ^ + ~ J)(^ j-i +ej)] + (z)A/tan0Q (E27) 

6c “ 2c 


Contributions from deflection terms 

Contributions from slope tenns 


(t=l,2,3,....,n) 


Equation (E27) is the degenerated form of equation (6c) for C,_j = C- = C 
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